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ABSTRACT 

This report describes an interpretive system which transforms the 650 
into a three-address, floating -decimal, general-purpose computer, 
primarily suited for scientific and engineering calculations. The system 
is complete in the sense that all mathematical, logical and input-output 
operations normally called for in such calculations can be performed 
within the system, i.e., without reference to the basic operation codes 
of the 650. The guiding principles in designing the system have been 
ease of use, as defined in the introduction, high speed of arithmetic and 
frequently used logical operations and full accuracy and range for the 
elementary transcendental functions. 

The report serves a dual purpose. It presents the external characteristics 
of the interpretive system to the potential user by means of detailed 
explanations accompanied by illustrative examples, assuming no previous 
familiarity with internally programmed machines. It also describes 
the internal structure of the system to the professional designer of such 
systems, enabling him to modify it to suit his particular needs or to 
borrow ideas or building blocks he may find useful. 

The program card deck for the system is available from the 650 Program 
Library. Requests should be addressed to: 

IBM 650 Program Librarian 
Applied Programming Publications 
International Business Machines Corporation 
590 Madison Avenue 
New York 22, New York 
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I INTRODUCTION 



I. 1. GENERAL DESIGN CONSIDERATIONS 



The use of most existing computing devices whose degree of automatic per- 
formance substantially exceeds that of a desk calculator entails certain problems 
not encountered in desk computing. To cope with these problems, one may incor- 
porate additional circuitry into the machine --this, indeed, appears to be the trend 
in recently announced commercially available machines--or, alternatively, one 
may program, in terms of the basic language of the machine, a system or 'super- 
language in terms of which the general user will program his problems. The user 
may consider the machine and the super -language as one entity, and no knowledge 
of the basic machine language is required of him. Before actual calculation, the 
programmer's instructions are translated by the machine into the basic language. 
If this translation or interpretation takes place each time an instruction is to be 
executed, rather than once for all at the beginning of a problem, the super-language 
is referred to as an interpretive language or system. Limitations in storage capac- 
ity may necessitate the choice of an interpretive system rather than a system of 
the once-for-all type in the case of most small or medium-sized computers. 

The designers of an interpretive system are faced with a very large number 
of decisions. To provide a basis of motivation for these decisions, it is convenient 
to list here, in somewhat arbitrary order, some of the above-mentioned problems 
which the present interpretive system proposes to solve. All of them may funda- 
mentally be measured in terms of total time spent by a programmer in learning to 
use the machine and in using it on a specific problem. In this sense, the "ease of 
use" referred to in the abstract above is implicitly defined by the list that follows. 
The price paid for the saving of programmer time is, of course, to be found in 
substantially reduced speed of operation. 

A. Scaling 

The storage medium- -paper- -normally used in desk computing places no 
practical restriction on the size of numbers or on the location of decimal points. 
In using a computer that automatically stores intermediate results in registers 
of fixed length and with the position of the decimal point fixed in advance, a great 
deal of time must in most cases be spent on estimating the range of all intermediate 
results to prevent errors due to overflow. The well-known way of avoiding this 
at the expense of a very substantial increase in the internal complexity of the arith- 
metic operations is to represent each non-zero number in floating -decimal form, 
i. e. , as a signed quantity whose absolute value lies in a fixed range, accompanied 
by an exponent of 10 or decimal point indicator. 

B. Length and Complexity of the Program 

Floating decimal arithmetic and frequently needed special functions could 



-4- 



be incorporated into a program written in the basic machine language in the form 
of a set of subroutines reached by a two-way transfer of control, ("calling se- 
quence" or "basic linkage") and there are indeed problems for which this is the 
best choice. In many cases—particularly in the case of relatively short prob- 
lems where the results are needed quickly- -a further reduction of the program- 
ming effort is desirable. This may be achieved by combining under single oper- 
ation codes, in an interpretive language, groups of steps in the basic language 
needed for performing frequently occurring tasks. For example, a single instruc- 
tion in which three locations are specified may be used for adding two arbitrarily 
located numbers and storing the result or a block of information of any length 
may be punched on cards as a result of a single instruction. In particular, the 
task of repeating a calculation a specified number of times, each time with appro- 
priate modifications, must be made easy, and the interpretive system described 
in this report goes as far as is believed possible in this direction by providing an 
order ("LOOP") with which simple cases of this task can be handled by a single 
instruction. 

To preserve the simplicity gained by introducing an interpretive system, the 
system must be made complete or self contained so that most problems can be 
conveniently programmed without reference to two different systems of operation 
codes, although, of course, leaving and re-entering the interpretive system should 
be made possible in order to provide the experienced programmer with complete 
flexibility. 

- C. Restrictions 

In desk computing, one cannot fail to notice if the argument for which a func- 
tion value is to be found in a table falls outside the range covered by the table or 
if in transcribing a set of numbers from one area on a piece of paper to another, 
overlapping area, some of the numbers to be transcribed are erased before tran- 
scription. A machine will avoid or detect such blunders only if programmed to do 
so, and as much as possible of this programming should be included in the inter- 
pretive system. For example, it is desirable that all mathematical functions in- 
cluded in the system be available for the full range of argument consistent with their 
definition and with the limitations imposed by the machine itself and that they are 
computed to the full accuracy of the number system used. Error stops indicating 
violations of unavoidable restrictions should be included to the fullest extent that 
space limitations permit. 

D. Program Testing 

The usefulness of a general -purpose computer or interpretive system depends 
decisively on the methods provided for testing ("debugging") new programs--for 
definitions and details, see Sec. VI. In the case of an interpretive system whose 
operating speed is only one order of magnitude above the speed of card punching, 
testing by means of a tracing routine included in the system compares favorably to 
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console testing, at least in the case of programmers whose familiarity with the 
machine is limited. Either of these methods is thought of as a tool normally used 
only when memory print-outs have been found insufficient. To facilitate testing 
by any of the methods mentioned and to keep programming as concrete as possible, 
the system described in this report assumes that the actual machine location of 
each instruction and stored number is assigned by the programmer. The system 
may, of course, be used in conjunction with regional or symbolic assembly programs. 

I. 2. CHANGES AND ADDITIONS 

Numerous minor changes suggest themselves when the system is viewed in 
the light of the experience gained in designing it; some of them are discussed in 
detail at the end of the report. Corrections of errors not yet revealed must be ex- 
pected. External changes and additions will undoubtedly be proposed after a period 
of use. The present system should, therefore, be considered primarily as a first 
version which each user may consider changing to better suit his needs. Comments 
and suggestions on both internal and external aspects of the system will be greatly 
appreciated. 
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II GENERAL INFORMATION 



II. 1. THE 650 

The IBM 650 is an electronic computer whose basic storage consists of 
a magnetic drum capable of holding 2000 words (numbers) of ten decimal digits 
and sign. The machine is internally programmed, i. e. , the program of instruc- 
tions which the machine isexpected to follow is kept on the drum, and the ma- 
chine automatically reads one instruction at a time from the drum, executes it, 
reads another instruction, and so on. Initially, the program is loaded onto the 
drum from punched cards but each instruction is loaded only once, although the 
machine may be expected to execute it many times in the course of a problem. 
Special orders are inserted into the program to cause repetition of prescribed 
sections the desired number of times. In many cases, the instructions exec- 
uted by the machine are changed or modified under control of the program between 
successive executions* This ability to modify its own instructions is one of the 
distinguishing characteristics of an internally programmed machine. 

In the basic model of the 650, all answers are punched by the machine into 
cards, which may be printed on a separate tabulator. 

The 650 is a general -purpose, fixed-decimal machine and any programmer 
may, with the aid of a detailed manual published by IBM, learn to use it as such. 
There are many large problems and many problems of a data processing nature 
for which fixed-decimal operation is definitely indicated, and the programmer is 
asked to give serious consideration to this alternative for all but the very smallest 
problems, since the gain in machine time over floating -decimal operation (ex- 
plained below) may be as high as 10:1. The machine -language programmer may 
relieve himself of many tasks by using the interpretive system for loading, punch- 
ing, calculation of special functions, etc. , provided 1000 storage locations suffice 
for his problem. (See TR OUT, Sec. IV. 1) 

This report describes a system which enables the programmer to use the 
650 as a floating -decimal machine, without being familiar with the fixed-decimal 
mode of operation. Beginning with the next section, all statements will concern 
the system rather than the 650 itself, but it should be borne in mind that if any- 
thing in the system appears restrictive from the viewpoint of a particular appli- 
cation, — storage capacity, speed, card form, word length, etc. --total or partial 
use of basic 650 coding may be the answer. 

II. 2. THE INTERPRETIVE SYSTEM: STORAGE; DATA AND INSTRUCTION FORM 

When the interpretive system is in use 999 ten-digit storage locations, 
numbered 001-999, are unrestrictedly available to the programmer for storing 
instructions and data. The location 000 has a special use ("previous result") 
which will be discussed below. 
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Throughout the system, numbers upon which mathematical operations are 
performed are stored and used in so-called (normalized) floating -decimal form, 
which will be defined as follows: The number zero is written as ten zeros with a 
plus sign ("machine zero"). Any number A other than zero is expressed as 



A=± A, 



io a i 



where 1^ Aj«<. 10 and -BO^ai^. 49. In the machine, A is written as the pair 
£(Aj, a), where a = a^-H 50 and Ai is an eight-digit number with seven decimal 
places. The "machine exponent" a is a two-digit (positive) integer located at the 
right end of the number. Non-zero numbers A not in the range lO'SO <^ | ^KlO^O 
cannot be correctly used in the machine, and some of the mathematical operations 
will give an error stop if the result would fall outside this range (see STOPS). Num- 
bers loaded into the machine must also be in the form prescribed above, unless 
special precautions are taken (see Sec. X. 2). 



The system instructions are signed ten-digit numbers of the following 



+ 


°1 


A or 2 


B 


C 



form: 



Here, 0j is a one-digit operation code and B and C are three-digit addresses. The 
three-digit quantity "A or 2 " is interpreted as an address A if 1 / and as an 
operation code O2 if 0^ ~ 0. The sign of the instruction is used in connection with 
the LOOP order (see LOOP OPERATIONS). The only difference between the mu- 
tually exclusive 0j and O2 operations is that all operations which require three 
addresses have been designated 0^, all others, O2. 

An example will illustrate how the addresses are used in a program: To add 
the number stored in register 200 to the number stored in register 201 and store 
the result in 500, the operation code 0^=1 (ADD) is used and the instruction reads: 
I 200 201 500. To take the square root of the number in 200 and store the result 
in 500, the 2 -operation 2 = 300 (SQRT) is used: 300 200 500. As will be shown 
in later sections, it is also possible to call out an instruction stored in memory 
and operate upon it, e.g. , increase one of the addresses in it. In storage, no dis- 
tinction is made between instructions and data so that the programmer is free to 
use any memory location for storing an instruction or a number as he sees fit. 

To facilitate explanations, the following notation will be used: The ten -digit 
quantity jvhose storage location has the address A — will be denoted by A; anal- 
ogously, B will denote the contents of location B. C denotes the result of a calcula- 
tion about to be stored in location C. 

In addition to being stored at C, the result C of any mathematical operation 
(i. e. , arithmetic operations or special functions) and of MOVE 000 and CONS 
(see READING) is automatically stored in the special location 000. If this result 
is needed on the next step, (or, more generally, before it has been replaced by 
the result of a subsequent mathematical operation or MOVE 000 or CONS) calling 
it out by using 000 as an A-address will reduce the execution time in the case of 
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the arithmetic operations. Also, time will be saved in any mathematical operation 
by using 000 as a C-address when C will be needed only on the next step. Execu- 
tion times are discussed in detail in a later section, but it should be emphasized 
that timing considerations only affect the running time of a problem, never the 
correctness of results. All locations are accessible at any time. 

Special addresses for obtaining frequently needed numbers, such as, and 
1, are not provided by the system. The programmer should load such numbers 
into locations of his own choosing. 
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Ill MATHEMATICAL OPERATIONS 



III. 1. ARITHMETIC OPERATIONS 

The operations will be introduced in an order chosen to facilitate the learn- 
ing process. Later, a concise summary of all operation codes will be given. 
Alphabetic operation codes are listed in addition to the numerical ones merely to 
facilitate programming; they are not introduced into the 650 and need not be used 
at' all. 



The result of each arithmetic operation is rounded. If the result is zero, 
a machine zero is given, i. e. , the machine exponent will be 00. An error stop 
occurs if the result of a multiplication or division would fall outside the range of 
the floating -decimal number representation; another error stop detects attempts 
to divide by zero (see STOPS). 



Numer. 
code 

0i = 1 



A list of the arithmetic operations follows: 

Alpha, 
code 



ADD 



Function 



Add Jin floating -decimal form) the num- 
ber A stored at A to the number B 
stored at B, store the result C at C 
and 000. Abbreviated: 



0i ss 2 



SUB 

MPY 
DIV 

NGMPY 
III. 2. SPECIAL FUNCTIONS 



1= 3 



X = 4 



°1= 5 



A + B = C 

Subtract: A - B = C 
Multiply: A . B = C 
Divide: A / B = C 
Multiply negatively: -A . B = C 



The system is intended to give eight-digit accuracy (i. e. , an error less than 
1 in the eighth digit) in computing the special functions included whenever the input 
makes this accuracy possible. For~trigonometric functions of an argument exceed- 
ing one revolution and for logarithms of numbers near 1, loss of accuracy follows 
from the mathematical properties of the respective functions and stops (which may 
be bypassed by the setting of a console switch) are provided when this loss exceeds 
two digits. For small values of the argument, where an eight-digit, fixed-decimal 
representation of the sine or arc tangent would contain leading zeros, the floating - 
decimal representation would normally introduce meaningless digits at the right 
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end. To reduce this nuisance to a tolerable level and also make possible trig- 
onometric calculations with extremely small arguments, the formulas sin x = x 
and arc tan x = x are used for lx\<, . 0025 and |x|< -. 001, respectively. Those 
interested will find the methods of computing the special functions described in 
Section XI. 3. 

Aside from the limitations imposed by the above mentioned inherent loss 
of accuracy and by the floating -decimal representation of the result, no restric- 
tions apply to the natural range of the argument for the special functions. Error 
stops will prevent attempts to take the square root of a negative number or the 
logarithm of a non-positive number. 

The special functions (or, more precisely, elementary transcendental 
functions) are: 



Numer. 
2 = 300 



2 = 301 



2 =: 302 



2 = 303 



0? = 304 



o = 305 



2 = 350 



2 = 351 



2 = 352 



2 = 353 
2 = 354 



2 = 355 



Alpha. 
SORT 
EXP E 
LOG E 
SIN R 
COS R 
ART R 
ABS 
EXP 10 
LOG 10 
SIN D 
COS D 
ART D 



Function 



B = C 
e B ~ C 



log B = C 
& e 



sin B = C, B in radians 
cos B 5= C, B in radians 
arc tan B = C, C in radians, ic|<CTT/2 
IbI = C 



10 B « C 



log 10 B s C 

sin B = C, B in degrees 

cos B = C, B in degrees 

arc tan B = C, C in degrees, |c|< 90 



Subdivisions of a degree are expressed decimally, not in minutes and sec- 
onds, i 

III. 3. MOVE 000 

In many cases (particularly in connection with the use of subroutines) it may 
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be convenient to be able to call out a number B from B and deposit it in C, as 
well as in 000, without the time-consuming use of a floating -decimal arithmetic 
operation. This is accomplished by the logical operation 0=9 ("MOVE") with 
A sr 000. The normal use of MOVE with A £ 000 is described in Sec. IV. 4. 

III. 4. AN EXAMPLE 

For the benefit of anyone with no previous computer experience, a simple 
example illustrating the use of the mathematical operations will be inserted here. 
Suppose that, as a part of a program which is assumed to be already on the drum, 
it is desired to evaluate the function 



f(x) = 



sin x 



/ 



1+ e 



-x? 



Here, x in radians is assumed to be in storage register 500 and the constant 1 in 
600. The quantity e~ x3 is to be stored in 501, and f(x) in 502. A program might 
look as follows: 

A or O2 B C Comments 

500 500 000 x 2 

000 500 000 -x 3 

3 
301 000 501 e~ x , store in 501 

.3 



Alpha. 


°1 


MPY 


3 


NGMPY 


5 


EXP E 





ADD 


1 


SQRT 





SIN R 





DIV 


4 



000 600 000 l+e"* x 

300 000 400 yjl + e" x3 , store 

temporarily 



303 500 000 sin x 

000 400 502 f( x ), store in 502 

The extensive use of the "previous result" address, 000, is worth notingo 



-12- 



IV LOGICAL OPERATIONS 



IV. 1. TRANSFER OPERATIONS 

Suppose the machine has been instructed (see LOADING) to begin a program 
by executing the instruction stored in, say, location 101. When this execution is 
completed, the machine will automatically execute instruction 102, then 103, etc. , 
until told by the program to do otherwise. Operations whose primary function is 
to influence either the order in which instructions are executed by the machine or 
the selection of stored data upon which the instructions make the machine operate 
will be called logical operations. A simple example of such an operation is 2 = 
203, "Transfer Control". If in the sequence 101, 102, 103 above, instruction 103 
should read "TR 203 000 080", the next instruction executed by the machine would 
be 080 instead of 104. This may be expressed by saying that "control was trans- 
ferred to 080". The B»address was ignored in this case. The transfer of control 
may be made to depend on the result of calculations (mathematical or logical) in 
which case a "conditional transfer" is said to occur. Logical operations --con- 
ditional or unconditional --are needed whenever several blocks of instructions, 
located on various parts of the drum, are to be tied together to form a program, 
whenever it is desired to repeat a calculation several times, etc. 

For simplicity in grouping, the following list of transfer operations includes 
two (UNC STOP and NOOP) whose transfer function is of a degenerate nature. In 
a first reading, it may be advantageous to omit the TR SUBR and TR OUT operations, 



Numer. 



Alpha. 



Function 



o 2 = 000 



UNC STOP 



Stop unconditionally. The machine 
stops regardless of the setting of 
console switches (see CONSOLE) 
and displays 9999 on the address 
lights and B on the display lights. 
This operation should be used only 
where it is intended to discontinue 
the execution of the program, since 
a continuation of the program cannot 
be effected by a simple depression 
of the PROGRAM START key (see 
STOPS). The C-address is ignored 
but should be filled, e.g. , with zeros 
(see LOADING). 



2 = 200 



COND STOP 



Stop conditionally and transfer. The 
machine stops if the PROGRAMMED 
STOP switch on the console is in 
the STOP position. The number 
1120 is displayed on the address 
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Numer. 



Alpha. 



Function 



lights and B on the display lights. 
When the PROGRAM START key 
is depressed, control is trans- 
ferred to C. If the PROGRAMMED 
STOP switch is in the RUN position, 
control is transferred to C with- 
out stopping. 



Caution: If the PROGRAMMED 
STOP switch is on RUN, the stops 
for loss of accuracy in sine, co- 
sine and logarithm and the stop in 
the CONS operation will not occur. 

This operation may be used for stopping at check points in the early running 
stages of a problem, with the option of avoiding the stops during later runs. 



0? = 201 



TR SGN 



2 sr 202 



TR EXP 



Transfer on sign. Control is trans- 
ferred to C if the result of the last 
mathematical operation or MOVE 
000 or CONS is negative, to B if 
it is non-negative, (i. e. , zero is 
regarded as having a plus sign). 

Transfer on exponent. The exponent, 
c, of the result of the last mathe- 
matical operation or MOVE 000 or 
CONS is compared to B (the leading 
digit of B should be 0). Control is 
transferred to C if c > B. If c< B, 
control proceeds to the next instruc- 
tion. 



This operation is particularly suited for the summation of series where terms 
are to be added until they have a prescribed number (50 - B) of leading zeros. For 
example, in order to return to instruction 080 only as long as the absolute value of 
the previous result is . 0001 or greater, one would write "TR EXP 202 046 080". 
This saves a time-consuming floating subtraction preceding the test. The TR EXP 
operation is also intended to take the place of the TR ZERO operation found in most 
systems. Due to the accumulation of small errors during a calculation, it is unwise 
in most cases to expect a result to be exactly zero to eight figures; here a TR EXP 
with a suitably chosen B may prevent a never-ending repetition of a part of a program. 



2 ss 203 



TR 



Transfer. Control is transferred to 
C, i. e. , the next instruction exe- 
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Numer. 



Alpha. 



2 = 204 



TR SUBR 



Function 

cutedwill be the one stored at C. 
The B address is ignored but 
should be filled, e. g. , with zeros. 

Transfer to subroutine. The C- 
address of the instruction located 
at C is set equal to B, whereupon 
control is transferred to C. The 
sign of the instruction at C is made 
positive. (For an elucidation and 
applications see SUBROUTINES. ) 



2 =■ 205 



TR OUT 



2 = 454 



NOOP 



Transfer out. Control is transferred 
to C and the instruction stored there 
is executed in the basic language of 
the machine (i. e. , outside the inter- 
pretive system). When an instruc- 
tion address 1095 is given in the 
basic language, control is returned 
to the interpretive system beginning 
at the instruction following the TR 
OUT. The B-address of the TR OUT 
instruction is ignored but should be 
filled, e.g. , with zeros. The pro- 
grammer in basic language must be 
careful not to use locations above 999, 
which are occupied by the interpretive 
system. 

No operation. Control proceeds to 
the next instruction. The B- and 
C-addresses are ignored but should 
be filled, e.g., with zeros. 



This operation is likely to occur in connection with tracing (see TRACING, 
particularly the ST TR ERAS operation) or when a superfluous instruction has 
been deleted from a program. 

IV. 2. LOOP OPERATIONS 

A highly repetitive character is required of any problem to be economically 
handled on an automatic computer. In certain instances, such as Newton's iteration 
procedure for the solution of equations, a repetitive process or "loop" is con- 
veniently programmed, merely using conditional transfer operations. In many cases, 
however, some of the instructions to be repeated must be slightly modified in a 
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systematic way before each new repetition. For example, in the evaluation of a 
linear expression £ a i ^ with the a£ and the Xi stored in blocks of consecutive 

i= 1 
locations, the addresses of a i and X£ must be increased by 1 each time a new term 
is to be computed. To facilitate programming of this kind, the system provides 
two methods of so-called address modification. The simpler—but less general-- 
of these methods employs a special counter called the loop box, which is stored 
in a location normally inaccessible to the programmer. If an instruction carries 
a minus sign, the current contents of the loop box will be added to the instruction 
(in fixed-decimal arithmetic and without regard to the sign) before it is executed. 
If, for example, the instruction - 1 531 600 901 is given and the loop box contains 
+ 009 000 009, the instruction actually executed by the machine would read 1 540 
600 910. The original instruction remains unchanged in its storage location. At 
the end of a calculation, an 0^ instruction called LOOP enables the programmer 
to increase the contents of the loop box by 1 in one or several address positions and 
to transfer control back to the beginning of the calculation. Hence, the calculation 
may be carried out repeatedly, each time with different addresses used in the exe- 
cution of instructions with minus signs. A test provision included in the LOOP order 
stops the repetition after a specified number of executions and resets the loop box 
to zero for future use. An example will be given after the following list of LOOP 
operations. 



Numer. 
2 = 100 



Alpha. 
LOOP A 



Function 

Loop on A. The contents of the 
loop box are increased by 001 
000 000. In other words, the A- 
segment of the loop box is in- 
creased by 1. After the increase, 
the A- segment of the loop box is 
compared to the B -address of the 
LOOP instruction. If the A- seg- 
ment is less than B, control is 
transferred to C. If the A- seg- 
ment is equal to B, {or greater, 
which will never be the case in 
normal use) the loop box is reset 
to zero and control proceeds to 
the next instruction. 



0, = 010 



LOOP B 



Loop on B. Analogous to LOOP A, 
with the B- segment of the loop 
box now being increased and com- 
pared to the B -address of the LOOP 
instruction. 



0, = 001 



LOOP C 



Loop on C. Analogous to LOOP 
A, with the C- segment of the 
loop box being increased and 
compared to B. 
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Nurner. 



2 = 110 



Alpha. 
LOOP AB 



o 2 = 101 



0? = Oil 



o 2 = 111 



LOOP AC 
LOOP BC 

LOOP ABC 



Function 

Loop on A and B. Analogous to 
LOOP A. The A- and B -segments 
of the loop box are increased by 
1 and the A-segment is compared 
to B. 

Loop on A and C. Analogous. 

Loop on B and C. Analogous. The 
B -segment is used for the com- 
parison. 

Loop on A, B and C. Analogous. 
The A-segment is used for the 
comparison. 



To illustrate the use of a LOOP order, consider the evaluation of the linear 
20 
expression L(x) = r^ a. x-, where the a. and the x. are stored in memory. In 

i =r 1 

choosing storage locations for numbers, it is wise to plan in advance how they are 

to be used in the program. In this case, since the a- and the x. are to be reached 

using the LOOP operation, it is advantageous to store them in blocks of consecutive 

locations, say the a i in 800 + i and the x i in 900 + i, (i s 1,2, . . . , 20). Suppose 

L(x) is to be stored in 700. For simplicity, assume that register 700 contains zero 

at the beginning of the calculation and that the loop box has been reset. The entire 

program for this calculation might be written as follows: 



Instr. No. 

101 
102 
103 
104 



Alpha. 



Sign. 







A or 0. 



B 



MPY — 3 

ADD + 1 

LOOP AB -+• 

Next instruction in the problem. 



801 


901 


000 


000 


700 


700 


110 


020 


101 



Note that the B -address of the LOOP order simply indicates the number of times 
the arithmetic calculation is to be performed, including the first time when the 
addresses are actually unmodified (modified by adding zero). The practice of 
starting the instruction numbering at, e.g. , 101, rather than 001 facilitates later 
additions to the beginning of a program. 

The loop box is automatically reset at the beginning of a new problem (see 
LOADING), and whenever a transfer out of a loop is effected by a loop order (as 
stated in the definitions above). Hence, the resetting of the loop box need not 
concern the programmer under normal conditions. If the need for resetting the 
loop box should arise, however, this is easily done by giving, e. g. , the order 
LOOP A with the B -address 000. According to the definition of LOOP A, this will 
cause control to proceed to the next instruction with a resetting of the loop box. 
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The C-address is irrelevant in this case. This situation would arise if control 
were transferred out of a loop in the middle of it by one of the conditional transfer 
operations. 

It is worth observing that a LOOP operation may be advantageously used in 
some cases where address modification is nof involved, simply to repeat a sequence 
of steps a prescribed number of times, e. g. , each time adding a fixed increment to 
a parameter. In such a case, any one of the loop orders could be chosen, (see 
EXECUTION TIMES, however) and no negative instructions would occur. 

The advantages of the loop-box method are its simplicity and high speed and 
the fact that the original instructions remain unchanged in memory. It is limited 
by the fact that there is only one loop box and hence, all instructions to be modified 
are modified in the same way. To handle situations more complicated than this, 
the system provides a set of operations described in the next section. 

IV. 3. ADDRESS CHANGE OPERATIONS 

Many problems can be completely programmed without the use of address 
change operations, and for someone approaching the field of internal programming 
for the first time, it might be advantageous to ignore this section until the need 
for more general logical operations arises. 

The functions of the address change operations are: (a) To increase or 
decrease a designated address of an instruction in storage by any given amount; 
(b) To set such an address to a given value (without reference to its previous value); 
and (c) To transfer control as a result of comparing such an address to a given num- 
ber. 

There are nine 2 -operations among the address c hang e_ ope rations. In each 
of these, the B-address gives the location of the instruction (B) to be changed and 
the C-address is the amount of change. For example, suppose the instruction 
600 750 005 (using the operation 2 - 600, ADD A) is given and suppose location 
750 contains the instruction I 320 400 000. Then the A-address, 320, of this in- 
struction will be increased by 005 and the resulting instruction 1 325 400 000 stored 
back in 750. Similarly, if 050 750 333 were given, (using 2 = 050, SET B) the 
instruction in 750 would be changed to read I 320 333 000. In brief: 



Numer. 



2 = 500 



Alpha. 
SET A 



Function 

Set the A-address. The A-address 
of the instruction (B) specified 
by Bis set equal to C. 



°2 = 



050 



SET B 



Set the B-address. The B-address 
of the instruction (B) specified by 
B is set equal to C. 



o 2 = 



005 



SET C 



Set the C-address. The C-address 
of the instruction (B) specified by 
B is set equal to C. 



-1 



Numer. 



2 = 600 



Alpha. 
ADD A 



Function 

Add to the A-address. The_A- 
address of the instruction (B) 
specified by B is increased by 
C. 



0. = 060 



ADD B 



Add to the B -address. TheJB- 
address of the instruction (B) 
specified by B is increased by 
C. 



0o = 006 



ADD C 



Add to the C-address. The C- 
address of the instruction (B) 
specified by B is increased by 
C. 



2 = 700 



SUB A 



o =r 070 



2 = 007 



SUB B 



SUB C 



Subtract from the A-address. 
The A-address of the instruction 
(B) specified by B is decreased 
by C. 

Subtract from the B -address. 
Analogous to SUB A. 

Subtract from the C-address. 
Analogous to SUB A. 



The sign of the instruction being modified remains unchanged and does not affect 
the outcome of the modification. Attempts to increase an address beyond 999 or 
decrease it below will result in erroneous operation not prevented by error stops. 

Three l -operations, TR A, TR B and TR C, complete the set of addresj, 
change operations. In each of them, the A-address specifies the instruction (A) 
to be called out and the B-address is the constant to which a specified address is 
to be compared. In case of inequality, control is transferred to C. For example, 
if the instruction 6 750 325 200 (using 0j = 6, TR A) is given, control will be trans- 
ferred to 200 if the instruction in 750 reads 1 320 400 000 but control will proceed 
ahead if 750 contains 1 325 400 000. Summarizing: 



0, = 6 



TR A 



Transfer on the A-address. The 
A-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 
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Numer, 



Alpha. 



Function 



0, = 7 TR B Transfer on the B-address. The 

B -address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 

l = 8 TR C Transfer on the C-address. The 

C-address of the instruction (A) 
specified by A is compared to B. 
Control is transferred to C if 
they are unequal but proceeds to 
the next instruction if they are 
equal. 

As an introductory example, the summation in the section on LOOP OPERA- 
TIONS will be programmed again using address change methods. This would be 
an inefficient choice in an actual problem, but it will best illustrate the difference, 
as well as the analogy between the two methods. It is again assumed that register 
700 contains zero at the start, but the steps analogous to the resetting of the loop 
box will be included. 



Inst. 



Alpha. 



Sign 



0. 



A or 0. 



B 



101 


SET A 


+ 





500 


103 


801 


102 


SET B 


-h 





050 


103 


901 


103 


MPY 


+ 


3 


[ ] 


[ ] 


000 


104 


ADD 




1 


000 


700 


700 


105 


ADD A 


+ 





600 


103 


001 


106 


ADD B 


+ 





060 


103 


001 


107 


TR A 


+ 


6 


103 


821 


103 


108 


Next ins true 


tion in the t 


)roblem. 









The brackets in the A- and B-addresses of instruction 103 are used to indicate that 
these addresses are variable and will be supplied by the program before the instruc- 
tion is executed, hence what is written there when the program is loaded into the 
machine is irrelevant. At the end of the program when instruction 108 is reached, 
memory location 103 will contain + 3 821 921 000. It is assumed that the summation 
just programmed is part of a larger problem in which it is used repeatedly. This - 
is the reason for the SET A and SET B instructions. If 801 and 901 were simply 
loaded into their respective positions in instruction 103 initially, the summation 
would be performed correctly the first time it is used, but the next time when the 
summation is called for, instruction 103 would read + 3 821 921 000 and erroneous 
calculations would result. The SET instructions could, of course, have been in- 
serted after the completion of the summation, restoring instruction 103 to its proper 
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value for the next application, but this procedure is not recommended because it 
makes it more difficult to restart the problem from the beginning without reload- 
ing the program in case of an interruption (e. g. , error stop) during the loop. 

A more realistic example of the use of address change methods would be a 
calculation involving more than one summation index or parameter. Then, one 
of the fast and convenient LOOP orders would normally be used in the "inner loop", 
i. e. , the loop occurring most frequently, with address change operations con- 
trolling the "outer loop" or loops. Suppose, for example, that it is desired to 
calculate S* ~ l ® a..x. for j = 1, 2, ... 5, where the a are stored in 800 + 10 + i 

i = l 
(i.e., a,, is in 811. a,, in 812, etc.; a^, in 821, a„„ in 822, and so forth), the 



ll is in 811, a 12 in 812, etc. ; a~, in 821, a in 82 
3 + i, and the S. are to be storedm 700 4- j. ft will 



x- in 900 + i, and the S~.~are to be stored in 700 4- j. It will be assumed that register 
500 contains zero. For completeness, the setting of all variable addresses to their 
initial values for repeated use of the summation program will be included. 



Instr. 


Alpha. 


Sign 


°I 


A or 0^ 
2 


B 


C 


Comments 


101 


SET A 


+• 





500 


104 


811 


Set variable 
' addresses to their 


102 


SET C 


+ 





005 


107 


701 

J 


initial values 


103 


MOVE 


+ 


9 


000 


500 


400 


Set register 400 
to zero 


104 


MPY 


mmm 


3 


Z 1 


901 


ooo"" 


"Inner loop" 


105 


ADD 


+ 


1 


000 


400 


400 


r i.e., summation 


106 


LOOP AB 


+ 





110 


010 


104 


on i 



107 



MOVE 



+ 



000 



400 



C 3 



Store the result 



108 


ADD A 


+ 





600 


104 


0101 


Increase addresses 


109 


ADD C 


+ 





006 


107 


OOlJ 


for next repetition 
in the outer loop 
(j-loop) 


110 


TR C 


+ 


8 


107 


706 


1031 


Test for end of 
j -loop 



111 Next instruction in the problem. 

A superficial examination of this program might suggest that only 1/5 of the 
program is devoted to actual arithmetic calculation ( I ), but it should be observed 
that in terms of the number of instructions executed by the machine when one 
complete summation is performed, the arithmetic ones are still in the majority, 
and in terms of execution time they comprise about 3/4 of the program. 
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In programming problems involving several loops, it may be helpful to 
consider the structure of a loop in terms of four phases: 

1. Initialization. Where addresses in the loop are set to their ini- 
tial values, registers used for summation are set to zero, etc. . 
The automatic resetting of the loop box and the fact that ad- 
dresses remain unchanged in memory tend to reduce the ini- 
tialization when the loop is controlled by a LOOP operation. 

In the summation problem above, steps 101 and 102 constitute 
the initialization for the outer loop, step 103 is the initialization 
for the inner loop. Notice that step 103 is repeated as a part of 
the outer loop. 

2. Execution. Comprising the mathematical operations of the loop, 
as well as any logical operations associated with a loop inside 
the one being executed. Above, the execution of the inner loop 
consists of steps 104 and 105 and the execution of the outer 
loop consists of 103-107. 

3. Modification. Where addresses, parameter values, etc. , are in- 
creased or decreased. The modification of the inner loop above 
is included in the LOOP instruction. The modification of the 
outer loop consists of steps 108 and 109. The position of the 
modification in the program in relation to the execution and test 
is frequently subject to choice. 

4. Test. Determining whether the loop is completed or further rep- 
etition is required. The LOOP instruction includes the test for 
the inner loop and step 110 is the test for the outer loop. 

Note: It is important to write loops in such a way that all ini- 
tializations are performed by the program, not by loading. If 
this rule is not followed, it will not be feasible to restart the pro- 
gram during testing or after a machine stop without reloading. 
For example, if a register is used for summation, it should be 
reset before being used in the summation loop by moving zero into 
it from another location, not by loading zero into it from a card. 

Many programmers find it helpful in programming a large problem to draw 
a block diagram or flow chart with one box representing each phase of each loop 
and arrows connecting the boxes showing the flow of control. 

The address change operations, particularly the SET operations, are fre- 
quently useful in non- repetitive situations as well. An example of this will be 
found in the section on SUBROUTINES. 

If a program appears to require a large amount of address modification 
and particularly, if this occurs because a quantity whose address is subject to 
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change is needed in many places in the execution of a loop, it may be advantageous 
to write the execution largely in terms of fixed addresses and perform the mod- 
ifications by moving data. Instruction 107, in the example above, illustrates this 
in a simple way: If the registers 700+ j themselves had been used in the summation 
process, (step 105) both the B- and C-addresses of step 105 would have required 
modification in the outer loop, as well as the C-address of instruction 103. For 
cases where several numbers are to be moved at the same time, a more general 
MOVE operation than the MOVE 000 used so far is available and will be described 
in_the next section. 

IV. 4. MOVE 

The MOVE operation is defined as follows: 



Numer. 
0j = 9 



Alpha. 
MOVE 



Function 

Move. If A ^ 000, the block of 
A consecutive words beginning at 
B is moved into the set of A con- 
secutive locations beginning at 
C. The words in the original 
locations are not destroyed, 
except where the two regions 
overlap. The number in location 
000 ("previous result") is not 
affected when C ^ 000. Both 
C y B and C ^ B are permissible. 
An error stop occurs if C 4- A — 
IS: 1000. If A= 000, the word 
(B) specified by the B -address is 
moved into location C and into 000, 
It also remains in location B. 



MOVE with A = 000 differs from MOVE with A= 001 only in that the execution time 
with A= 000 is shorter and that the previous result location is affected. 

Note: If a number is to be moved from location B into 000 for use in a TR 
SGN or TR EXP operation on the next step, MOVE 9 001 B 000 must not be used, 
since these transfer operations work strictly according to their definitions (see 
Sec. IV. 1. ). The correct instruction would be MOVE 9 000 B 000. (Internally, 
these transfer operations inspect a duplicate "previous result" location rather than 
000!) 

As an example, suppose xj is in 701, x 2 in 702, . . . , x 5 in 705 and the instruc 
tion MOVE 9 005 7 01 703 is given. Then x x will be found in 703, x 2 in 704, . . . , 
x^ in 7 07, after execution. 

In conclusion, it should be pointed out that the use of the logical operations 
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is by no means restricted to the straightforward functions for which they are 
primarily intended. The programmer will find innumerable ways of increasing 
the efficiency and elegance of his programs by unusual applications, particularly 
of the address change operations. As a weird example, suppose it is desired 
to multiply the numbers located in registers 1, 4, 9, 16, 25, 36, . . . , 400 (! ) 
by a constant located in 600 and store the results in 501, 502, 503, . . . , 520: 



instr. 


Alpha. 


Sign 




°1 


A or 0^ 


B 


C 


898 


SET A 


4- 







500 


900 


001 


899 


SET C 


+ 







005 


901 


003 


900 


MPY 


— 




3 


[ ] 


600 


501 


901 


ADD A 


+ 







600 


900 


C ] 


902 


ADD C 


+ 







006 


901 


002 


903 


LOOP C 


+ 







001 


020 


900 


904 


Next instruction in 


the ] 


problem. 









-24- 



V INPUT-OUTPUT OPERATIONS 



V. 1. CARD FORM 

By a card form is meant a specific assignment of card columns to form fields 
for data, instructions, identification, etc. , in connection with a given program or 
interpretive system. In the 650, information is transmitted to and from cards 
through a control panel, and anyone whose needs call for a special card form can 
adapt it for use in connection with the interpretive system merely by simple control 
panel wiring. For most needs, the following card form, associated with the inter- 
pretive system, is likely to be found adequate. At this point, only brief definitions 
of the card fields will be given for reference in subsequent sections where their 
use will be explained in detail: 



Columns 



Definition 



1-4 
5-6 
7-9 

10 

11 
12-21 

22 
23-32 

33 
34-43 

44 
45-54 

55 
56-65 

66 
67-76 
77-79 

80 



Card number 

Deck number 

Location 

"Word count 

Sign of word 1 

Word 1 

Sign of word 2 

Word 2 

Sign of word 3 

Word 3 

Sign of word 4 

Word 4 

Sign of word 5 

Word 5 

Sign of word 6 

Word 6 

Problem number 

Tracing identification 



The same card form is used in all input-output operations, as well as in trac- 
ing. Both instructions and data are signed ten-digit words and are entirely indistin- 
guishable in connection with input-output operations. 

V. 2. PUNCHING 

At any point in the problem, the machine may be ordered to punch into cards 
the contents of any set of memory locations, together with appropriate identification. 
In some problems, it may be desirable to punch out answers one at a time, perhaps 
together with the values of relevant parameters; in others it may be preferable to 
punch out a large amount of information at less frequent intervals. There are also 
cases where it is advantageous to punch out instructions r In connection with testing 
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ioW^^ TESTING) in ° rder to «"»i»e a program interrupted at a chosen 

a deck of ' cardTllW 1^ ^ ng ' (See L ° ADING > in - der to "duce the size of 
a deck of cards. All of these ends are served by the following instruction: 



2 =: 410 



PCH 



Punch cards. The block of consecutive 
words beginning at B arid ending at C 
(inclusive) is punched into cards. Five 
words and a word count of 5 are punched 
into each card but the last, whose word 
count will be the remainder when C-B + 
1 is divided by 5. On each card, the 
location from which word 1 was punched 
is punched into columns 7-9. The 
words in storage are not destroyed. A 
cumulative count of the number of cards 
punched during the problem (i. e. , since 
LOADING) is punched into columns 1-4. 
The problem number (see LOADING) is 
punched into columns 77-79 and zero is 
punched into columns 6 and 80. An error 
stop occurs if B > C. 

™„-£ " ^a K^ ir ! d . t0 PUnCh SlX W ° rdS t0 a Card ' thiS ma ^ be done b y adding a 
special card behind the punching deck (see LOADING). This card should have an 

TrZl 00 Z^Znn - 1% Y n C ° 1UmnS 6 ' 9 ' l ln C ° 1Umn 10 > a 12 P— h in column 
11, and 00 0006 0000 m column 12-21. 

The punched cards are likely to be used for one (or both) of two purposes (in 
addition to possible processing on other equipment): The information on them may 
be printed on a tabulator or they may be loaded (or READ) into the 650 at a later 
time. Details of the printing will not be given here, since they depend on the char- 
acteristics of the tabulator, but the printing form may be assumed to be rou g hlv 
identical with the card form with proper spacing between words. (Suggestions on 
tabulator wiring are given in Sec. XI. 5. ) It is assumed that the suppression of the 
superfluous words punched into the last card, if its word count is not 5, will be 

i 6 A°cTu d °,V he tabulat ° r contro1 P anel - « ^is is not feasible, it may be done in 
the 650 by adding three cards to the punching deck. For details, see Sec. XI. 1. 

Selective spacing between lines in printing may be accomplished in several 
ways, even though no operation in the 650 is provided for this purpose A brief 
discussion will be given here, since spacing considerations may affect the use of 
the PCH operation in programming. Through the setting of switches on the tab- 
ulator, a choice of any of the following spacing alternatives may be provided: 

(a) Single or double spacing. 
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(b) Spacing between every n lines (with n chosen by wiring, nor- 
mally, e.g., nslO). 

(c) Spacing after any line whose word count is less than the word 
count of the preceding line. 

(d) Spacing before any line whose location number has a units digit 
smaller than the units digit of the location number of the pre- 
ceding line. 

Alternative (c) is suited for the printing of information punched from fairly 
large blocks of locations by one PCH order. Spacing will occur after each block, 
unless the block length is a multiple of 5, which can be avoided by programming. 
Alternative (d) is intended for information punched repeatedly from the same set 
of locations and provides the option of spacing when the loop is interrupted, e. g. , 
for changing a parameter value, 

V. 3. LOADING 

When a program has been written, and careful inspection reveals no further 
errors, it is key punched into cards following the card form given in Sec. V. 1. To 
reduce to a minimum the number of errors to be found with the aid of the 650, the 
cards should be run through a verifier operated by another person or, alternatively, 
key punched independently by two operators and compared on a reproducer. The 
programmer has the option of specifying the number of words to be punched to a 
card: Punching 5 or 6 to a card will keep the program deck small from the outset 
and eliminate the need for condensing the deck on the 650 later. Punching one word 
to a card is felt by some programmers to facilitate changes. Each card must have 
in columns 7-9 the location into which word 1 is to be loaded, and in column 10 the 
number of word3 to be loaded from the card into consecutive locations. Columns 
1-6 and 77-80 are not read by the 650 (except that the problem number is read from 
the last card, see below) and may be used by the programmer as he deems best. 
Each column of each field to be used by the machine must contain one and only one 
punch and an error stop is provided to enforce this rule. A 12-punch is used for 
plus, an 11- or x -punch for minus and a 0-punch--not a blank column- -for zero. If 
the word count is less than 6, unused word fields and sign columns may be left 
blank. No distinction is made between data and instructions in key punching and load- 
ing. 

LOADING is the process of feeding data and instructions into the machine at 
the beginning of a problem. If the previous user of the 650 was not using the inter- 
pretive system or if there is any reason to doubt that the system is correctly 
stored on the drum, the program deck should be preceded in loading by a deck which 
loads the interpretive system (in 51. 9 seconds) into the memory locations above 999. 
Before the program deck, the programmer may also place a Reset Memory Card, 
which will (in 6. 3 seconds) reset each of the memory locations 001-999 to minus 
zero. (This is useful in connection with the punching out of sections of memory in 
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testing. ) Immediately behind the program deck--no blank cards are used in the 
card reader in connection with this interpretive system- -the programmer places 
one of two nine -card decks to inform the machine whether he wants normal opera- 
tion or TRACING described in a later section. (If he knows that he wants the 
same mode of operation as the previous user, he can omit these cards but the 
gain is only 2. 7 seconds. ) Last, he must place a so-called transfer card with a 
zero punched in column 10, the problem number in 77-79 and the location of the 
instruction at which the program begins in columns 7-9. The word fields on this 
card may be left blank. 

The loading program automatically resets the loop box, the card counter 
(see PUNCHING) and location 000 to zero. 

The order in which the program cards are loaded is irrelevant, unless the 
same location is loaded into from more than one card, in which case the last such 
card, of course, determines the contents of the location. This may occur in 
connection with changes of a temporary nature, which may be placed at the end of 
the deck and later removed, leaving the program in its original form. In the deck 
which loads the interpretive system, the order of the cards must be preserved, 
and an error stop is provided to insure this, thereby ascertaining that no part of 
the system is missing. 

In summary, a complete deck to be loaded must contain: 

System deck (173 cards) 

Reset Memory card (optional) 

Program deck 

Mode -of- ope ration deck (9 cards) 

Transfer card 

The control console of the 650 need be of almost no concern to the user of 
the interpretive system under normal conditions. He must only make sure that 
all switches on the console are set in a fixed manner required by the system, and 
these settings will now be listed without any description of the function of the 
switches. Certain ways of using the console are described in the sections on READ 
ING and PROGRAM TESTING. 

Switches Settings 

Storage entry 70 1951 1333 + 

Programmed stop Stop (see COND STOP) 

Half cycle Run 

Address selection 1338 (see STOPS) 

Control Run 

Display Upper Accumulator 

Overflow Stop 

Error Stop 

To start a problem, the deck to be loaded is placed in the card reader, and 
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the following keys are depressed in order: 

(1) COMPUTER RESET 

(2) PROGRAM START 

(3) START 



(on the console) 
(on the console) 
(on the card reader ] 



When the last card leaves the hopper, the machine stops and the key labelled 

(4) END OF FILE (on the card reader) 

is depressed. If the deck has been correctly put together, the execution of the pro- 
gram will then start automatically. 

The program deck may be run out at any time after loading by depressing the 
START key, unless a READ instruction is contained in the program. Blank cards 
should be inserted into the PUNCH hopper and the START key on the punch side 
depressed. 

To make the 650 produce a condensed program deck in case the program was 
originally key punched one instruction to a card, a PCH instruction should be given 
at the very beginning of the program. This instruction may be bypassed during 
subsequent executions of the program merely by changing the location number on 
the transfer card. 



V.4. READING 

In some problems, particularly in applications of a data processing nature, 
it may be desirable to read information into the machine during the execution of 
the program without manual interference. This is accomplished by the READ 
operation: 



2 = 400 



READ 



Read cards. The block of consec- 
utive storage locations beginning 
at B and ending at C (inclusive) 
is read into from cards. The ad- 
dress B must appear in the location 
field on the first card, as well as 
in the READ instruction, and the 
location field on each card follow- 
ing must contain the sum of the 
word count and location on the 
previous card. The sum of the 
word counts of all cards to be 
read must be C-B+ 1. Violations 
of these requirements, which have 
been included for the programmer's 
protection, will result in error 
stops. 
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The cards to be read should be placed in the hopper of the card reader imme- 
diately following the transfer card (no blank cards). 

The decisions made with the aid of conditional transfers and other logical 
operations are normally based on criteria predetermined by the programmer and 
incorporated into the program. If the programmer wishes to influence the pro- 
gram during its execution, e. g. , on the basis of a result displayed on the console 
in connection with a COND STOP instruction, he may do so using the CONS opera- 
tion: 

2 — 401 CONS Read console. The machine stops 

if the PROGRAMMED STOP switch 
is on STOP. Zero is displayed on 
the display lights and 1131 on the 
address lights. When the PROGRAM 
START key is depressed, the num- 
ber entered on the STORAGE 
ENTRY SWITCHES is stored in 
location C and in 000 (the "previous 
result' r location). If the PRO- 
GRAMMED STOP switch is on 
RUN, the storing takes place with- 
out a stop preceding. The B -ad- 
dress is ignored but should be 
filled, e.g. , with, zeros. 

An example of an application of CONS might be the feeding in of an "educated 
guess" for a starting value in connection with the solution of algebraic equations. 
Another application, involving only the storage entry SIGN switch, might be to 
continue a program until another user is ready to take the machine, at which time 
a change in the SIGN switch setting, interpreted by a TRSGN operation, causes the 
program to punch out intermediate results for later restart. 
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VI PROGRAM TESTING 



VI. 1. MEMORY PRINT-OUT 

The choice of methods for testing ("debugging") a program by comparing 
results of machine calculation to known quantities or to results of independent 
calculations by other means is governed by the relative availability of machine 
time and programmer time. If machine time is freely available, testing with the 
aid of the control console is highly efficient, as well as instructive and enjoyable, 
as soon as a certain facility for operating the console has been acquired. Par- 
ticularly in the case of small problems, the method of tracing --where a card is 
punched for each instruction executed, showing all numerical and logical quantities 
associated with the execution- -may be the most desirable in that it gives an almost 
certain clue to the difficulty within a predictable, if not very short, period of 
machine time and allows the programmer to study the material at his leisure. 
The method most economical of machine time and yet frequently sufficiently 
illuminating is that of memory print-out. It might be suggested that on most prob- 
blems in a busy but not heavily over-loaded installation, the methods be used in 
the order reverse to that in which they were mentioned here. Some directions 
for their use will now be given. 

The memory print-out method simply consists of inserting temporarily into 
the program at one or several suitably chosen points PCH orders (see PUNCHING) 
calling, for the punching of blocks of information- -data or instructions --which, when 
printed on the tabulator, will give a picture of the progress of the program. Since 
1000 words may be punched 5 to a card in two minutes, it is not out of the question 
to punch out the contents of every register used in a problem- -including all the 
instructions --several times. To get the most benefit from this method, the pro- 
grammer should, in any problem that does not threaten to fill the entire available 
memory, avoid using the same storage location for storing different quantities at 
different times whenever feasible, so that as many partial results as possible are 
preserved for the memory print-outs. Whenever a test case of a problem is run, 
even if memory print-out is not chosen as the primary testing method, it would 
certainly be advisable to make the last instruction of the test deck punch out the 
entire memory used. A flexible alternative would be to have scattered through the 
program CONS — TR SGN combinations which transfer control to a PCH order if 
the storage entry sign switch is turned to minus. 

Temporary instructions may be inserted into a program in two ways: Either 
they areincluded in the normal sequence of instructions when the program is 
initially written and replaced either by NOOP instructions (see Sec. IV. 1) or by 
transfer to the next non-temporary instruction when no longer needed,or else one 
of the regular instructions of the program is replaced by a TR to a vacant location 
L, the regular instruction is placed in L,, the temporary ones inL+1, L+ 2, etc. , 
and at the end of this temporary sequence a TR back to the normal program is given. 
In either case, the temporary instructions may (as suggested in LOADING) be kept 
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as a separate deck at the end of the program deck, eliminating any changes in the 
main program deck and simplifying bookkeeping. 

VI. 2. TRACING 

If the tracing deck of nine cards is loaded with the program deck, (see LOAD- 
ING) the machine will automatically start tracing from the beginning of the program, 
as specified by the transfer card. Before the execution of each instruction, a card 
with the following information will be punched: 

Columns Definition 

1-4 Card number (cumulative) 

6 Zero 

7-9 Location of the instruction 

about to be executed. 
10 Six 

11-21 The instruction as stored in 

memory. 
22-32 The instruction as modified 

for execution (i. e. , with the 
loop box added if minus). 
33-43 The contents of the loop box. 

44-54 A if A #000, zero if A « 000. 

55-65 B 

66-76 The contents of location 000 

(i. e. , the result of the last 
mathematical MOVE 000 or 
CONS operation). 
77-79 Problem, number. 

80 Eight (used by the tabulator 

for automatic selection of a 
different printing form for 
trace cards). 

The punching rate will be 100 cards per minute except in the case of very time-con- 
suming operations, such as, the moving of a large block of information. The ad- 
vantage of punching the trace card before execution is that information will be 
punched for an instruction whose execution is interrupted by an error stop. In the 
case of instructions (such as LOOP or TR EXP) whose B -address does not refer 
to a memory location, the quantity B is irrelevant. Tabulator wiring to suppress 
the printing of B in such cases can be provided if sufficient selection equipment is 
available. The PCH operation is bypassed when the machine is operating in the 
tracing mode, i. e. , PCH is equivalent to NOOP. 

If a program is too long to be traced in its entirety or if this is unnecessary, 
selective tracing may be effected by using the following operations: 
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2 = 450 



START TR 



CL = 451 
2 



STOP TR 



Start tracing. If the nine -card 
tracing deck has been loaded, 
the "machine will start tracing 
from the next instruction. If it 
is already tracing, it will con- 
tinue to trace. The B and C 
addresses are ignored. If the 
deck for normal operation has 
been loaded, START TR will 
be equivalent to NOOP. 

Stop tracing. If the machine is 
tracing, it will discontinue trac- 
ing immediately. If it is not 
tracing, STOP TR will be equiv- 
alent to NOOP. The B and C 
addresses are ignored. 



452 



ST TR ERAS 



Start tracing and erase itself. 
If the tracing deck has been 
loaded, the machine will start 
tracing from the next instruc- 
tion. If it is already tracing, 
it will continue to trace. If the 
deck for normal operation has 
been loaded, . tracing will not 
begin. In all cases, the ST TR 
ERAS instruction will be replaced 
in memory by a NOOP (02= 454) 
during its first (and only! ) execu- 
tion. The B and C addresses are 
ignored. 

The bypassing of the PCH operation is in effect as long as the trace program is on 
the drum and is not affected by the selective tracing orders. To make PCH opera- 
tive, the nine-card deck for normal operation must be loaded. 

The purpose of the ST TR ERAS operation is to make it possible to trace the 
repetitive steps of a loop either once or twice and then stop tracing until the loop 
is completed. To get the steps traced once, one may place the pair STOP TR, ST 
TR ERAS at the beginning of the repeated portion of the loop; to get them traced 
twice, one places this pair of instructions at the end immediately preceding the 
test. As a specific example, suppose it is required to trace twice the steps of the 
loop programmed in the section on LOOP OPERATIONS and suppose vacant loca- 
tions are available from 900 up. Assume that the machine is tracing as it enters 
the loop. The original program reads as follows: 
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801 


901 


000 


000 


700 


700 


110 


020 


101 



101 MPY — 3 

102 ADD -+■ i 

1 03 LOOP AB + 

1 04 Next instruction in the problem. 

The following instructions could be added as a temporary deck at the end of the pro- 
gram deck: 

102 TR + o 203 000 9 00 

900 ADD + i 000 7 00 700 

901 STOP TR + 451 000 000 

902 ST TR ERAS + 452 000 000 

903 TR + 203 000 103 

Notice that the TR instruction gets loaded into 102 after the regular program, replac- 
ing the ADD instruction, as explained in LOADING. This example is, of course 
unrealistic in that selective tracing would hardly be needed for testing such a simple 

VI. 3. CONSOLE TESTING 

Testing with the aid of the control console requires some familiarity with the 
internal structure of the interpretive system (see Sec. XI) and with the basic language 

T IT ^^lTalV e l ting iS m ° re attractive on the 65 than on most machines due 
to the ADDRESS STOP feature: If the CONTROL switch is turned to the ADDRESS 
STOP position, the execution of the program will proceed at electronic speed until 
the address set up on the ADDRESS SELECTION switches is reached. At that point, 
the machine stops, and the contents of various registers may be examined on the 
display lights or the program may be continued manually one step at a time Alter- 
natively, the program may be punched out on cards at this point by merely feeding 
m one card with a PCH instruction, going into any vacant location, followed by a 
transfer card specifying this location. Console testing, in connection with the inter- 
pretive system, is likely to be needed only in exceptional cases. 

The ADDRESS STOP feature of the 650 may be used in conjunction with a spe- 
cial address stop transfer card when it is desired to start tracing from a certain 
instruction N in the middle of a program after running at full speed up to that point. 
Uhis may, of course, alternatively be accomplished using the tracing operations 
described m Sec. VI. 2, but then the value of N must be decided upon in advance 
and the proper program changes key punched. ) The procedure is as follows- Set 
the ADDRESS SELECTION switches to N and turn the CONTROL switch to ADDRESS 
STOP. Load as usual and run until the machine stops at the instruction N. (For 
details on possible earlier stops see below. ) Then set the CONTROL switch to 
RUN and load the tracing deck followed by the address stop transfer card. Tracing 
will begin immediately and the first instruction traced will be N. 

In choosing N it must be remembered that the loop box and location 000 are 
reset to zero when the tracing deck is loaded. If this restriction is inconvenient, 
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it can be circumvented by placing a special card in front of the tracing deck. The 
card counter and the problem number are also reset to zero, unless the tracing 
deck has been modified to prevent it. 

If the CONTROL switch is kept in the ADDRESS STOP position when the pro- 
gram deck is loaded, one stop will occur when location N is reset by the memory 
reset card and another when the programmer's instruction is loaded into N. Also, 
stops may occur before instruction N is reached in the program, if N is referred 
to in an ADDRESS CHANGE or MOVE operation (but not if N is one address in a 
conditional transfer instruction and control is transferred to the other address). 
After each stop, operation will resume when the PROGRAM START key is depressed. 
If the CONTROL switch is left in the ADDRESS STOP position during tracing, two 
stops will occur each time N is referred to (and one if N is the B -address of a 
transfer instruction). 

The program can be continued at full speed (punching mode) after a period of 
tracing by following the procedure described above with the punching deck in place 
of the tracing deck. 

The address stop transfer card has 69 1976 1952 24 1061 1098 in columns 
1-20 and a 12-punch in each of columns 1, 10 and 20. The special card for by- 
passing the resetting steps in loading has 69 1953 1952 24 1278 1953 70 1951 1344 
in columns 1-30 and a 12-punch in each of columns 1, 10, 20 and 30. (See Deck 7, 
Sec. XI. 7. ) 

If the value of N has been decided upon in time to get it key punched into a 
regular transfer card, (Sec. V. 3) this card may, of course, be used in place of 
the address stop transfer card in the procedure described above. 

A programmer familiar with the internal structure of the interpretive system 
will be able to think of many other cases where special needs can be met using 
maching language cards ( ,r load cards"). 
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VII SUMMARY OF OPERATION CODES 



Oj OPERATIONS 



2 OPERATIONS 







Page 






Page 






Page 


Num. 


Alpha. 


Ref. 


Num. 


Alpha. 


Ref. 


Num. 


Alpha. 


Ref. 





GO to 


8 


000 


UNC STOP 


13 


300 


SQRT 










200 


COND STOP 13 


301 


EXP E 




1 


ADD 


10 


201 


TR SGN 


14 


302 


LOG E 




2 


SUB 


10 


202 


TR EXP 


14 


303 


SIN R 




3 


MPY 


10 


203 


TR 


14 


304 


COS R 




4 


DIV 


10 


204 


TR SUBR 


15 


305 


ART R 




5 


NGMPY 


10 


205 


TR OUT 


15 


350 
351 


ABS 
EXP 10 




6 


TR A 


19 


100 


LOOP A 


16 


352 


LOG 10 




7 


TR B 


20 


010 


LOOP B 


16 


353 


SIN D 




8 


TR C 


20 


001 


LOOP C 


16 


354 


COS D 










110 


LOOP AB 


17 


355 


ART D 




9 


MOVE 


23 


101 


LOOP AC 


17 














011 


LOOP BC 


17 


400 


READ 


29 








111 


LOOP ABC 


17 


401 
410 


CONS 
PCH 


30 
26 








500 


SET A 


18 














050 


SET B 


18 


450 


START TR 


33 








005 


SET C 


18 


451 


STOP TR 


33 








600 


ADD A 


19 


452 


ST TR ERAS 


33 








060 


ADD B 


19 














006 


ADD C 


19 


454 


NOOP 


15 








700 


SUB A 


19 














070 


SUB B 


19 


- 












007 


SUB C. 


19 
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VIII STOPS 

Error circuits in the 650 will stop the machine on attempts to use invalid 
information, such as, that represented by blank columns or double punches, as 
well as on several kinds of machine malfunctioning, and will indicate on the control 
console the nature of the error. If this occurs during the loading of a new deck, 
the cards should be examined. In other cases, a note should be made of the indica- 
tions on the console, and the procedure that led to the stop should, if possible, be 
repeated exactly in order to determine whether the error is systematic in nature. 

All stops, which are part of the interpretive system, will now be listed. 
Conditional stops will occur only if the PROGRAMMED STOP switch is set to 
STOP. On a conditional stop, the PROGRAM LIGHT in the OPERATING section 
of the console will be on and no lights in the CHECKING section should be on. The 
program will continue if the PROGRAM START key is depressed. On an uncon- 
ditional stop, the STORAGE SELECTION light in the CHECKING section will be on. 
Normally, operation should be discontinued after an unconditional stop and changes 
made in the program in order to avoid the stop. Alternatively, the program may 
be continued by having a transfer card (see LOADING) in the card reader, specifying 
the instruction to which control should proceed when the COMPUTER RESET and 
PROGRAM START keys are depressed. 

The location of the interpretive system instruction xxx on which the machine 
has stopped, may be determined by displaying the contents of location 1098 on the 
console. The display lights will show 60 Oxxx 1107. This process, called "monitor- 
ing !t , may be performed as described in the 650 manual or, alternatively, by setting 
the storage entry switches to 60 Oxxx 8000 and depressing the COMPUTER RESET, 
PROGRAM START and PROGRAM STOP keys. 

If, in an exceptional case, it would be advisable to proceed to the next instruc- 
tion after an unexpected unconditional stop, this may be done manually as follows: 

(1) Set the CONTROL switch to MANUAL. 

(2) Check that the ADDRESS SELECTION switches are set to 1338. 

(3) Depress the COMPUTER RESET key. 

(4) Depress the TRANSFER key. 

(5) Set the CONTROL switch to RUN. 

(6) Depress the PROGRAM START key. 

As a result of this procedure, zero will be stored at C and 000 before the next instruc- 
tion is executed. If this is not desired, the ADDRESS SELECTION switches should 
be set to 1095 in step (2). To repeat the same instruction (on which the stop occurred) 
the switches are set to 1098. 

There is an alternative manual procedure for restarting after an unconditional 
stop which is simpler in the case of frequent use but is not recommended in general 
because it requires changing the setting of the STORAGE ENTRY switches. They 
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are used in LOADING and must be set back to their normal positions for the next 
user: 

(1) Set the STORAGE ENTRY SWITCHES to 00 1951 1338+- (or 
00 1951 1095 + if zero is not to be stored or 00 1951 1098 + 
to repeat). 

(2) Depress the COMPUTER RESET key. 

(3) Depress the PROGRAM START key. 
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CONDITIONAL STOPS 



Address Lights 
1120 
1131 
1715 

1835 



■Normal Cause 

Programmed COND STOP. (Display lights show B. ) 

CONS (Check STORAGE ENTRY switch setting. ) 

Loss of accuracy in SIN (Exponent of B exceeds 52) 
or COS. 

Loss of two or more digits of accuracy in LOG. 



UNCONDITIONAL STOPS 



2222 



r 



i 



MOVE with C + A - 1 >_ 1 000. 
PCH with B > C + 1. 



.READ with incorrect loc. or word count. 



3333 
4444 



DIV with B = 0. 
SQRT with B < 0. 



5555 



MPY with result out of range. 
DIV with result out of range. 



6666 



^EXP with result out of range. 

LOG with B <. 0. 

SIN with exp. of B exceeding 58. 
yCOS with exp. of B exceeding 58. 



7777 



Cards missing or out of order in the system deck 
being loaded. 



9999 



Programmed UNC STOP (Display lights show B 
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IX EXECUTION TIMES 

The execution times listed in this section are based on the standard 650 drum 
speed of 12, 500 r. p. m. They represent approximate theoretical estimates derived, 
in the case of the mathematical operations, from simple assumptions regarding the 
distribution of the numbers to be operated upon. For example, the part A, of a 
floating -decimal number A = Aj . 1 0^1 is assumed to be uniformly distributed 
between 1 and 10, although in physical problems there are reasons that favor a 
logarithmic distribution; extremely small and extremely large exponents are con- 
sidered very unlikely, etc. It is further assumed that the programmer has chosen 
storage locations on the drum without regard to timing, ignoring the fact that in 
the case of some operations the execution time will be a few milliseconds shorter 
for numbers stored in certain sections of memory. Some, but not nearly all, of 
the time estimates have been verified by tests. 

It should be stressed that the estimates of execution times are needed only 
for making comparisons or estimates of running time for problems or for choosing 
efficient ways of programming and will never affect the result of an operation. In 
comparing these estimates to estimates given for other interpretive systems or 
subroutines, it is important to verify by sample calculations or machine tests that 
the assumptions are realistic. 

To minimize the size of the table, the execution times listed refer to a basic 
case and corrections to be added in other cases are given at the beginning of the 
table. 



-40- 



650 INTERPRETIVE SYSTEM 
ESTIMATED AVERAGE EXECUTION TIMES IN MILLISECONDS 



(a) If A^ 000, add 7. 2 ms. for ADD and SUB, 6. 3 ms. for MPY, NGMPY and 
DIV*. 

(b) If C £ 000, add 6. 1 ms. for all mathematical operations, MOVE 000 and 
CONS*. 

(c) If the instruction has a minus sign, add 4. 8 ms. for all operations. 

(d) If, after a TR EXP or LOOP operation, control will proceed to the next 
instruction rather than to C, add 4. 8 ms. 



ADD 
SUB 

MPY 

DIV 

NGMPY 



TR A 
TR B 
TR C 



MOVE 



65. 7 


UNC STOP 


28. 8 


65. 7 


COND STOP 


29.8 


67. 2 


TRSGN 


19.2 


74. 3 


TR EXP 


24. 


67. 2 


TR 


19.2 




TR SUBR 


44.4 




TR OUT 


26. 


37. 3 


LOOP A 


24. 


37, 3 


LOOP B 


28.8 


42. 1 


LOOP C 


24. 




LOOP AB 


24. 




LOOP AC 


24. 


37. 7 


LOOP BC 


28.8 


40. 8 + 12A 


LOOP ABC 


24. 


(A = no. of 






words. ) 


SET A 


55. 3 




SET B 


55. 3 




SET C 


55. 3 




ADD A 


44.9 




ADD B 


44. 9 




ADD C 


44. 9 




SUB A 


44. 9 




SUB B 


44.9 




SUB C 


44.9 



SQRT 
EXP E 
LOG E 
SIN R 
COS R 
ART R 
ABS 
EXP 10 
LOG 10 
SIN D 
COS D 
ART D 

READ 



CONS 
PCH 



206 
197 
202 
192 
187 
238 

33.2 
187 
207 
240 
235 
271 

One card: 

101+ 14n 

(n= no. of words. ) 

Succeeding cards: 

300 each. 

28.8 

One card: 

163+ 12. 5n 

(n = no. of words. ) 

Succeeding cards: 

600 each. 



START TR 


28. 8 


STOP TR 


24.0 


ST TR ERAS 


38.9 


NOOP 


24.0 


TRACING 


600 per card, 


LOADING 


300 per card 



*(See next page for footnote. ) 
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"'Those who are particularly interested in time considerations may wish to know 
the exact increments on which the weighted averages in (a) and (b) are based: 

In ADD and SUB, 4. 8 ms. if 17 < A < 41 (mod 50) 

9.6 ms. if 1 < A< 16 or 42 < A< 49 (mod 50) 

In MPY, NGMPY and DIV, 4. 8 ms. if A > 17 or A = 1 (mod 50) 

9.6 ms. if 2 <A< 16 (mod 50) 

In all mathematical operations, MOVE 000 and CONS: 

4. 8 ms. if 7 < C < 42 (mod 50) 

9. 6 ms. if 1 < G < 6 or 43 < C < 49 (mod 50) 

An easily remembered programming rule could be extracted from this in- 
formation: If locations between 17 and 41 (mod 50) are used for storing numbers, 
the increments given in (a) and (b) may be replaced by 4. 8 ms. 
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X SPECIAL TOPICS IN PROGRAMMING 



X. 1. SUBROUTINES; TRANSLATION 

A subroutine is a program expected to be of use as a part of the program in 
several problems or in several sections of the same problem. The mathematical 
operations in the interpretive system are indeed subroutines written in the basic 
language of the machine and reached through their operation codes, and anyone 
desirous of preparing an additional subroutine of this type may avail himself of a 
vacant operation code (see Sec. XI) and write the program in a part of the memory 
below 1000. 

Subroutines written wholly or partly in the interpretive language may be 
reached conveniently using the TR SUBR operation defined in Sec. IV. 1. Suppose 
the subroutine begins at 900 and ends at 9 35. Instruction 900 should read: SET 
C 005 935 C Hand instruction 935 should read: TR 203 000C 1 Now suppose 
the subroutine is needed at step 700 in a program, and when it has been used, 
control is to be transferred to 680. Instruction 700 should read: TR SUBR 204 
680 900. The TR SUBR operation will take the quantity ("return address") 680, 
place it in the C-address of instruction 900 and then transfer control to 900. Instruc- 
tion 900, in turn, places 680 in the C-address of 935, and when instruction 935 is 
reached at the end of the subroutine, it transfers control to 680 as originally desired. 
Hence, the programmer using the subroutine only needs to know the identifying 
number 9_0 ; the transfer of control to and from the subroutine is handled by the 
TR SUBR in conjunction with the two instructions 900 and 935 provided by the writer 
of the subroutine. Subroutines needing only one input number and giving only one 
result (such as, the evaluation of one Bessel function for a given value of the argu- 
ment) will normally assume the input to be in 000 and will deliver the result there; 
in the case of several numbers, specified locations normally within the block 
occupied by the subroutine would be used for input and/or results. Subroutines may, 
of course, be used inside other subroutines without restriction. 

If the locations occupied by a subroutine are needed for another purpose, e. g. , 
another subroutine in the same problem, the subroutine may be translated to a 
different set of locations by a translating program developed by Miss D. C. Leagus. 
When the subroutine is written entirely in the interpretive system, the programmer 
is required only to separate data and constants from instructions, and the translating 
program will automatically decide which addresses of each instruction are subject 
to translation. Machine language instructions may also be used in a subroutine to 
be translated, provided certain conditions specified by the translating program are 
adhered to. 

Subroutines for the solution of cubic equations and of systems of linear 
equations have been written at the Laboratories. 
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X. 2. UNNORMALIZED INPUT; TRANSITION BETWEEN FLOATING- AND FIXED- 
DECIMAL FORM 

Nearly all of the mathematical operations in the system assume that the float- 
ing decimal numbers to be operated upon are in the normalized form defined in 
Sec. II. 2, i. e. , that the leading digit is different from zero unless the entire num- 
ber is zero. In processing empirical data, key punching is often facilitated by 
permitting leading zeros and reproducing a constant exponent. Such unnormalized 
data may be used in the interpretive system provided the first operation in which 
it is used is ADD or SUB with operand exponents differing by less than 10. 

A special case of unnormalized input is that of a zero with a non-zero machine 
exponent. If such a zero is added to a non-zero number with a smaller exponent, 
a number of digits equal to the difference between the exponents are lost. Con- 
sequently, zero should be equipped with exponent 00 unless the programmer knows 
in detail how the zero will be used in his program. Special provisions in the MPY 
and DIV routines make it possible to use a zero with machine exponent 00 in tnem 
without danger of exceeding the exponent range negatively. 

The converse problem of producing unnormalized output, e. g. , for the print- 
ing of tables in fixed-decimal form or for calculations in machine language is easily 
solved at the expense of one digit. Suppose for example that the numbers N. to be 
"unnormalized" or "fixed" are known to be less than 10 4 and output in the form 
XXXX. xxx is desired. Add the number 1000000054 (i. e. , 10, 000. ooo) to N. if 
N^ > 0, subtract it from N. if N. < (using TR SGN) and punch. The output of the 
form -+- (10, 000 4- |Nj) , is ready to be printed on the tabulator with the leading 
1 and the constant exponent 54 suppressed by hammerlocks or wiring. If the num- 
bers are to be used in machine language, the 1 and 54 are shifted out. Rounding 
to a smaller number of digits is obtained by choosing the exponent of the additive 
constant (1000000054) correspondingly larger. 

X. 3. EXAMPLES 

In conclusion, two problems will be programmed in order to illustrate the use 
of many of the operations and methods described. 

First, suppose it is desired to evaluate the "error function", 
for a set of values x s a, a+& , a+2A, . . . , a+10A, using the RAND approximation 

(2> <^*cx)= I -Ca|n+a E n z +a 3 Tf+a^n 4 -t-a 5 n sr )<^'cx) ) 
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where 



(3) 



ri~l/Cl + px) ; £ p js a numerical constant). 



and 



(4) Q.'CX)=-£- 



e-x 



and to punch out the results as well as to store them for later use. 
of the polynomial inYt will be faster if (2) is written in the form. 



The evaulation 



(5) ^*Cx)=i-Cn(a 1 +TiCa 2 -t-n(a 3 +n(a 4 +na 5 )))))^ , (x> 



«*, 



To make it possible to use the LOOP order in evaluating Gf(x) this way, the 
coefficients a^ will be stored in consecutive locations in decreasing order. The 
LOOP program will be given a form applicable to an arbitrary polynomial by 
including a "dummy" coefficient a Q = 0. Storage locations will be chosen as 
follows: 



Location 



Contents 



101-119 

200 
201 
202 
203 
204 

221 
222 
223 

301 
302 
303 
304 
305 
306 

401-410 



instructions 




QfC*) J 



(cards 1-4) 



constants 
(card 5) 



temporary 

storage 

("erasable") 



coefficients 
in <£*C*) 
(card 6) 



results 
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The program might be written as follows; 
Alpha. Sign (^ A or B 



SET C 



MOVE 



+ 



4- 







005 



000 



114 



202 



C Comments 

401 Set address of. 

first <£*C*) 
221 First x is x = a 



103 


NGMPY 


4- 


5 


000 


104 


EXP E 


+ 





301 


105 


MPY 


4- 


3 


000 


106 


MPY 


4 


3 


204 


107 


ADD 


4- 


1 


000 


108 


DIV 


4- 


4 


201 


109 


MOVE 


4- 


9 


000 


110 


MPY 


+ 


3 


000 


ill 


ADD 


— 


1 


000 


112 


LOOP B 


4- 





010 


113 


NGMPY 


4- 


5 


000 


114 


ADD 


4- 


1 


000 


115 


ADD 


+ 


1 


221 


ii6 


ADD C 


4 





006 


117 


TR C 


4 


8 


114 


118 


PCH 


4 





410 


119 


COND STOP 


4 





200 



000 


000 


-xZ 


000 
200 


000 
222 




221 


000 


px 


201 


000 


1+ px 


000 


223 


n=i/ci+px) 


301 


000 


a 5 into 000 for LOOP 


223 


000 


prev. res.,Tl 


302 


000 


add next coeff* 


005 


110 


loop in polynomial eval 


222 


000 


polyn. '<p'C>Q 


201 


C 1 


Q (x)= 14 prev. res. 


203 


221 


x + A » next x 


114 


001 


next<£*(x) address 


412 


103 


test for end 


401 


410 


punch two cards 


221 


500 


end; stop, display last 
x, go to 500 on PRO- 
GRAM START 



500 Next instruction in the problem. 

An important remark should be made: If there is no shortage of storage 
locations and if the programmer does not mind writing a somewhat larger number 
of instructions, the running time for many problems can be decreased and the 
logic simplified by "unwinding' 1 the innermost loop, i.e., by writing out the 
mathematical instructions in the loop in a straight sequence instead of using the 
LOOP operation. In the present problem, a sequence containing five MPY and four 
ADD instructions could replace the instructions 109-112 and also eliminate the use 
of the dummy coefficient a . The execution time for the polynomial loop would be 
reduced by nearly 1/3 and the LOOP operation could be used to replace the ad- 
dress change operations in the outer loop. The polynomial evaluation accounts for 
about 1/2 of the total running time of this problem. In many large problems, the 
innermost loop consumes an even larger fraction of the running time, making it 
important to program the innermost loop efficiently even at the expense of apparent 
inefficiencies elsewhere. 



x 



V 



The second illustrative problem reads as follows: For a given set of numbers 
v =1, 2, . . . , 50, not necessarily equally spaced, the values of the Chebyshev 



polynomials T n (x v ), n«l, 2, . . . , 10, are to be computed using the recursion formulae 
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(6) T n+l Cx v )=dx v T n Cx v )-T Tl _,Cx v ) ; 

(T Q (x v ) si, T^ (x v )= x v ) and punched out in a compact form. 
In addition, the sum 

(7) 5 s ° CT,o<:Xv)3 2 ( X -xy) 
V=IYI-X V * 

(Xci = 1) is to be punched out and the operator is to be given the option of also 
calling for the punching of partial sums of (7) at any time. 

Storage locations will be assigned as follows: 

Location Contents 



050 
051 
052 



0) 
1) 



constants (card 7) 



095-120 instructions (cards 1-6) 



199 
200 
201 
202 

210 

300 

301 

302 

400 +■ v 
451 



The sum (7) and its partial sums) 
T (x v ) a 1 
T x (x v ) 5 x v 
T 2 (x y ) 



output 



Tiq( x v) 

^Xy 

[ T io< x J] 2 

[Tio(*v)] 2 /Y 

x v 
1 



"^ 



> 



J 

} 



temporary- 
storage 
{''erasable 1 ') 



input 

(cards 8 - 18) 



In addition, locations 1-8 will be used in connection with a trick in program- 
ming the LOOP. 

The program may be written in many ways. The following is not necessarily 
the best: 
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13 
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CM 
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LD 
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ft w 
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CM 
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CM 
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CM 
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o 
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rH rf 
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A number of remarks are called for, many of them of general applicability: 

(a) The C-address of instruction 101 will, during execution, run through 
the values 000-008, but the result of the instruction is always called out from 000 
on step 102. This trick makes it possible to use the LOOP BC operation instead 
of address change, which is normally required if different sets of addresses are to 
be modified during a loop. 

(b) The instruction numbering was arrived at by starting the preparation 
of the program at instruction 101 with the intention of later adding an unknown 
number of initialization steps preceding it. This speaks in favor of not starting 
a program at 001. 

(c) The stop which would normally occur each time the CONS instruction 
is reached may be bypassed when found superfluous without any sacrifice by 
turning the PROGRAMMED STOP switch to RUN, since no COND STOP, SIN or 
LOG operations (the only other ones involving a conditional stop) are used. The 
operator decision regarding punching of partial sums is made using only the sign 
switch of the STORAGE ENTRY switches. This switch does not influence LOAD- 
ING. 

(d) The quantity x is used so frequently that it was more economical to 
MOVE it into a fixed location than to apply address modification. The converse 
applies to x v -|- 1? which is used only once. 

(e) The constant 1 appears in three locations merely in order to simplify 
bookkeeping and loading, as well as changing the number of points x v in a later 



run. 



(f) An invaluable aid in determining whether the results of a calculation 
are correct is a mathematical identity which they must satisfy, and the programming 
of such checks is strongly recommended whenever it is possible. In the present 
problem, the identity 

1 

(iii) \ LiiJ dx = X 



V 



1 — x 



2 



is closely connected with the computation of (7) if the x v are distributed over the 
interval (-1, 1). 

(g) An alternative method of programming the outer loop, which would 
eliminate the address change operations at the expense of somewhat increased 
card preparation, would be to key punch the x v one to a card and give a READ 
order entering one x v at a time into a fixed location during the execution of the 
program. The difficulty arising from the need for x v + \ on step 109 is not in- 
surmountable. 
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XI INTERNAL STRUCTURE OF THE SYSTEM 



XI. 1. DETAILED DESIGN CONSIDERATIONS 

An expert examining the program at the end of this report will ask a number 
of questions about apparent duplication, about tight optimization in one routine in 
contrast to a lack of it in another, about the choice of operations and of methods 
of implementing them, etc. This section will attempt to answer some of these 
questions and also suggest a number of changes and additions that could be con- 
sidered for a second version of the system. Additional questions and suggestions 
from readers will be genuinely appreciated. 

In the early stages of system design, the following requirements were among 
those agreed upon, in addition to the general principles discussed in Sec. I. 2: 

(a) The arithmetic operations and those logical operations most likely to 
occur in inner loops (LOOP and certain TRANSFER operations) must be as fast 
as we know how to make them, regardless of the expense in storage. 

(b) The system must occupy at most 1000 memory locations. 

(c) The special functions must have full accuracy and unlimited range and 
most of them should be as fast as these requirements and available storage permit. 

(d) Optimum programming, (see the 650 Manual of June, 1955) in addition 
to being necessary for the attainment of (a) and (c), should be used locally in any 
program where the gain is significant but not at the expense of extensive rewriting 
of previously completed programs. 

(e) The programs must be so written that if the machine stops on any pro- 
gram step in a subroutine and control is transferred elsewhere before restarting, 
the subroutine, where the stop occurred, is left in a condition which assures correct 
operation the next time that subroutine is used. This implies that if a subroutine 

is used in more than one program, it must be initialized by each program rather 
than having a normal form used in one program and temporarily being changed at 
the beginning of other programs when needed there and then restored to normal at 
the end. 

(f) To facilitate changes, the individual programs (or "decks", 1-20, see 
Sec. XI. 7) that make up the system should be as independent of one another as they 
can be without excessive waste of storage. This requirement was not fully adhered 
to near the end of the programming task. 

As a result of these requirements and of some oversights in programming, 
there are a number of storage registers which could be made available without any 
loss in system performance and a number which could be freed at some sacrifice. 
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A brief guide for finding such registers will now be given followed by a number of 
suggestions for their possible use in a revised version of the system. 

The 6 vacant 2 ~code locations and the 11 vacant registers listed in deck 5 
are, of course, available. The only distinction between them is one of mnemonics 
in connection with the choice of operation codes. In addition, it appears possible 
to salvage 22 registers essentially without loss by the following substitutions, but 
a careful check followed by machine testing is advisable: 



)eck 


Card 


Loc. 


Replace by 


20 


113 


1801 


1848 


16 


48 


1240 


1138 


18 


103 


1896 


1138 


6 


70 


1360 


1160 


16 


6 


1230 


1338 


18 


77 


1887 


1137 


2 


32 


1058 


1955 


12 


60 


1639 


1289 


16 


30 


1244 


1245 


19 


103 


1702 


1103 


17 


63 


1480 


1980 


10 


36 


1331 


1358 


8 


32 


1166 


1241 


12 


56 


142* 


1674 


17 


22 


1485 


1285 


18 


97 


1842 


1504 


17 


29 


1495 


1297 


5 


11 


1252 


1952 


5 


12 


1255 


1955 


5 


13 


1260 


1960 


5 


31 


1277 


1977 


5 


23 


1283 


1983 



It is, of course, necessary to determine, by sorting on instruction and data 
address, all places where the locations listed are referred to. 

Registers that may be freed at a price in speed include, above all, nearly 
40 extra registers used in the arithematic routines in calling out A and B, splitting 
them up and storing the parts. This is done separately in each of decks 12, 13 
and 14 to accommodate minor differences that facilitate optimization. To combine 
these steps without any loss of time is a task which, if possible, would require 
re -optimization of a substantial part of the system. At the expense_of one revolution, 
they may be combined easily. Similarly, making the dissection of B common to 
all 2 - routines would result in a substantial saving at the expense of lost time in 
cases (such as LOOP and TREXP) where B is irrelevant. To make this dissection 
common only to those routines where it is needed would be less profitable. 
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At some sacrifice in external characteristics, registers may, of course, be 
freed in any number of ways. If, in tracing, the modified instruction (redundant 
but convenient) is omitted, seven steps are eliminated. The MOVE operation for 
A£ 000 is easily programmed in terms of LOOP BC and MOVE 000 and could be 
omitted, as could the special functions in degrees and to base 10 (or radians and 
base e, respectively). 



A number of suggestions for changes and additional operations will now be 
listed. Suggestions (1) - (3) use only the vacant registers and operation codes 
listed in deck 5 and can consequently be added to the system without difficulty at 
the option of any installation or individual programmer. For temporary use, they 
may be punched on separate cards and loaded after the system deck, in the case 
of (1) and (2) and after the punching deck, in the case of (3). Such cards should 
have an x-punch in column 5 and the four-digit location in columns 6-9. 



(1) Add an 0? -operation defined as follows: 



2 = 453 



SWITCH 



Transfer on switch. Control is 
transferred to C if the Storage 
Entry Sign Switch is set to minus, 
to B if it is set to plus. 



This operation bypasses the stop that would occur if the same function were 
programmed by a CONS and a TR SGN order. It might be particularly useful in 
connection with tracing when it is desired to start tracing after a certain amount 
of running time has elapsed or for following the progress of a calculation by 
occasional punching of intermediate results at the discretion of the operator. 



The coding for SWITCH consists of the instruction: 
1453 10 8000 1015 



Read console. Go to TR SGN 
routine. 



The execution time is 19. 2 ms. 

(2) Add an ? -operation called COUNT having the same counting and testing 
properties as the LOOP orders but using a counter independent of the loop box and 
not capable of modifying instructions. Its function can be duplicated, e.g., by a 
SET A, an ADD A and a TR A instruction. Its advantage lies in its speed and 
simplicity. The execution time is 24. ms. when control is transferred to C and 
33. 6 the last time when control proceeds ahead and the counter is reset. A formal 
definition follows: 



2 = 800 



COUNT 



The number standing in the counter 
is increased by 1. Its new value 
is compared to B. If B is greater, 
control is transferred to C. Other- 
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wise, the counter is reset to 
zero and control proceeds to 
the next instruction. The 
counter is also reset in loading. 

If COUNT is used extensively, an expansion of the tracing program to punch 
out the contents of the counter, e. g. , in place of the problem number, would seem 
desirable. 

The coding for COUNT reads as follows: 



Loc. 



Op. 



1800 


10 


1314 


10 


1323 


11 


1381 


46 


1337 


10 


1396 


21 


1391 


16 


1066 


20 


1356 


00 


1317 


00 


1194 


20 


1378 


24 



Data 

1356 
1317 
8002 
1337 
8001 
1356 
8002 
1356 

0000 
1000 

0000 
1356 



Instr. 




0000 
0000 

1378*1 
1178J 



Remarks 

Call out and increase 

the contents, N, of the counter. 

Test N+ 1- B. 

On — , store N •+• 1 in the counter; go to C 

(in TR SGN program). 

OnH- , reset the counter, go to General 

Interpretation. 

The Counter 
Constant 

Change in LOADING to reset 
the counter. 



Note: If the COUNT program is loaded separately, the card loading 
zero into 1356 (step 9 in the program) must be included. 

(3) Include in the punching program (deck 9) a routine that prevents unwanted 
numbers from being punched out when the word count is less than the normal max- 
imum. This can be done on a tabulator with sufficient selector capacity (see Sec. 
XI. 5). In the 650, it requires five locations and increases the execution time of 
the PCH order by 24. ms.when the word count is less than the normal maximum. 
The program, which may be punched on three cards, reads as follows: 



1949 


44 


1306 


1095 


1306 


20 


1980 


1307 


1307 


20 


1981 


1308 


1308 


20 


1982 


1309 


1309 


20 


1983 


1044 


1044 


20 


1984 


1973 



(4) Make room for the tracing program to be on the drum in parallel with 
the punch program, replacing the mode-of-operation deck (see LOADING) by an 
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x-punch on the transfer card or a setting of the storage entry sign switch. An ex- 
pansion of the loading program (about 7 steps) or of general interpretation would 
be needed, and the present overlap between tracing and punching is 34 registers. 

(5) If MOVE is omitted, except for A = 000, make this an O^-operation and 
use the vacant Oj-code for NGDIV. Alternatively, add a fast -operation, "NEG", 
identical with MOVE 000 except that it changes the sign of B. If NEG were available, 
however, it might be used in cases where, by slight reprogramming, a better pro- 
gram using NGMPY could be written. 

(6) Increase the number of logical operations, adding to the flexibility of the 
system and to the confusion of the beginner: Have a register called the !r address 
counter", addressed, e.g. , by 000 or by special operations and SET instructions 
referring to the address counter (as in 701 Speedcoding) where the present SET 
instructions refer to their own C-address. Have a set of TR A, TR B, TR C orders 
which automatically increase the address referred to by 1. These would have to be 
alternative to the address transfer orders in present use unless vacant 0, -codes are 
produced. 

(7) Make use of addresses now ignored in some operations. For example, 
in CONS, use the B -address to call out a number B for console display when the 
machine stops. In TR OUT, make B a "return address" similar to that in TR SUBR. 
In START TR, or a new tracing order supplementing it, let B (or C) designate the 
number of steps to be traced before an automatic discontinuation of tracing. 

(8) Make Program Loading reset the registers below 1000 to zeros, unless 
told not to by an x-punch on the first card being loaded. 

(9) Have a conditional stop, or an operation effecting such a stop, on loss 

of accuracy in ADD and SUB, analogous to those in SIN and LOG. In many problems, 
particularly in connection with tests, such loss is legitimate, however,and a stop 
undesirable. 

(10) Replace or supplement the present error stops by the punching of an 
"error card". 

(11) Introduce an operation similar to ST TR ERAS, perhaps replacing it, 
which will cause the machine to trace the first, second and last repetition in a 
loop. 

(12) Add another LD-STD pair (at no loss in time) to General Interpretation 
(see cards 26 and 27) making ADD and SUB, as well as MPY and DIV available as 
internal subroutines. 

(13) Cut the execution time of several subroutines, such as, the arc tangent 
program, by making minor rearrangements, usually involving the expenditure of 

a few additional registers. 
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(14) Add an 0~-operation, SPACE, which causes an x-punch to appear on the 
next card punched. 

(15) Interchange the functions of registers 1002 and 1702, causing the machine 
to stop sooner if a programmer accidentally attempts to continue upward from 
instruction 999. 

(16) Investigate whether a carry can ever occur on card 78, deck 18. If not, 
put the registers used on cards 78-82 to better use. 

(17) Replace or supplement the arc tan operation by an ] -operation, ARG, 
which gives the argument (angle) of the point whose coordinates are (A, B). 

XI. 2. RELATED SYSTEMS 

Several systems supplementing the present one suggest themselves: (a) A 
system of symbolic or regional programming where the machine assigns absolute 
addresses in connection with loading; (b) A system externally identical with the 
present one, or very nearly so, operating on complex numbers, probably with real 
and imaginary parts in 8-2 floating -decimal form; (c) A system externally identical 
with the present one, or very nearly so, operating on double precision floating - 
decimal numbers, e.g. , 16-4; (d) A system of "formula translation" or "automatic 
coding" (such as, the IBM Fortran for the 704) putting on the machine as much as 
possible of the burden of translation from a set of mathematical formulas to a 
program. 

XI. 3. NUMERICAL METHODS 

The study of numerical methods for calculation of the special functions included 
in the system was not nearly as exhaustive as would have been desirable and no 
claim to an optimal choice is made. 

The square root is computed by Newton's iteration method, 

where 1 \ B, £ 10, using the initial approximation 

(2) x Q s 1 4- .22 B 1 . 

The evaluation of the trigonometric and exponential functions is based on RAND 
approximations (see Approximations for Digital Computers by Cecil Hastings, Jr. , 
Princeton University Press) to sin ^ X and 1 X for 0^x< 1. Resembling the 
approximations obtainable by expansion in terms of orthogonal polynomials, these 
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approximations are in general somewhat more efficient than partial sums of Taylor 
series for a prescribed interval and accuracy, but it is not obvious that a further 
reduction of the argument followed by the use of a Taylor expansion could not have 
been better in the present case. For small x, as stated previously, the formula sin 
x a x (in radians) is used in order to retain significant figures. 

The logarithm and, arc tangent are evaluated from fixed-length partial sums 
of power series after preliminary reductions of the argument, since eight-digit 
RAND approximations were not available. For the logarithm of B,, 1< B, <^ 10, 
the substitutions 



(3) 



U= jk ; V=£ , for B,< e ; 



(4) U "L6? >V=U5 ; WB^e, 

(5) + s. ttrl 
x u+l ' 

are followed by the evaluation of 

The constants, VlTand e** °5 were arbitrarily chosen within the intervals that would 
lead to a minimal number of terms in (6). For x near 1, the logarithm is inherently 
less accurate than x since 

(7) d log X= ^ ^ ix 

and log x s= whereas x^ 1. No substitution comparable to sin x = x can alleviate 
this difficulty. The use of a second (^-operation for log (1 4- B) was considered 
but was rejected due to space limitations. This second logarithm could not be used 
to replace the present one for all values of the argument, since if the logarithm of 
a small number, say 10~ 10 is desired, the substitution, 

(«> I 4- B = icf ,0 

forced upon the programmer, yields B= -1 exactly (in the eight-digit system used) 
with all digits of the input lost. 



-56- 



_ For the arc tangent, the reduction (after the argument is restricted to 0<C 

B\l by the obvious properties of the function) is based on the formula 

(9) arc tan x= arc tan y +■ arc tan */ 

l + Xy 

which is merely the addition theorem for the tangent rewritten. With y= .6, the 
use of (9) gives the desired accuracy in 

(10) arc tan z - Z |~J _ **+ Z *_ , , ,_ Z!?l 

u 3 5 I! J 

with z =s -£—< for x^.28 and z = x for x <^.28. There is again some leeway in the 

choice of these constants. For small x, the substitution arc tan x = x is used to 
preserve significant digits. 

XI. 4. CONTROL PANEL WIRING FOR THE READ-PUNCH UNIT 

The control panel for the 533 Read-Punch Unit associated with the 650 is 
wired as follows: 

Col. 1, 1st Reading, to LOAD. 
R -H Sign, jackplugged. 
P -f- Sign, jackplugged. 

Col. 5, 1st Reading, to Pilot Sel. 1 X PU. 

Rd. Hold to PS1 Hold. 

Read Card C, Col. 6, to PS1 T. 

Read Impulse to PS1 N. 

PS1 C to Storage Entry C, Word 1, pos. 3 (from the left). 

Read Card C, Col. to Storage Entry C. 

7-9 Wd. 1, pos. 4-6 

10 Wd. 2, pos. 6 

11 Wd. 3, Sign 
12-21 Wd. 3, pos. 1-10 
22 Wd. 4, Sign 



67-76 Wd. 8, pos. 1-10 

77-79 Wd. 9, pos. 4-6 
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Read Impulse 12 Wds. 1, 2, 9, 10, Sign 

Wd. 1, pos. 7-10 
Read Impulse -^ Wd. 2, pos. 7-10 

Wd. 9, pos. 7-10 



Word Size Emitter to 



10 
8 



Word Size Entry C 

Wds. 3-8, 10 
Wd. 1 



7 Wd. 9 

5 Wd. 2 



Storage Exit C to 



Punch Card C, Col, 



Wd. 


1, 


pos. 


3-6 


Wd. 


2, 


pos. 


6 


Wd. 


3, 


Sign 




Wd. 


3, 


pos. 


1-10 


Wd. 


4, 


Sign 




Wd. 


8, 


pos. 


1-10 


Wd. 


9, 


pos. 


4-6 


Wd. 


2, 


pos. 


10 



Wd. 10, pos. 3-6 1-4 

6-9 
10 
11 
12-21 

67-76 
77-79 
80 

Double Punch and Blank Column Detection as available and desired. 
XI. 5. CONTROL PANEL WIRING FOR THE TABULATOR 

The IBM accounting machine or tabulator used for printing from the cards 
associated with the interpretive system may be expected to perform some or all 
of the following tasks: 

(1) Automatic selection of different printing forms (i. e. , zero control and 
spacing between items on a line) for data cards and trace cards. 

(2) Selective spacing between lines. 

(3) Suppression of unwanted words from cards with word count less than 
5 (or 6). 

Since there are many tabulator models, it is not feasible to provide a detailed 
wiring diagram in this report. Instead, suggestions of general applicability will be 
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given. 

Exact selector requirements depend on the characteristics of each machine. 
As an example, requirements on a 416 will be given: The printing of signs requires 
6 single -position selectors with X-pickup and 6 positions of 11-12 separation (either 
special attachments or 6 positions of a selector transferred by an 11-1/2 impulse). 
Task (1) requires a digit selector (which may be put to duplicate use in task (3)) and 
34 selector positions with digit pickup (delayed pickup). Task (2) requires only 
3 one-position selectors on the 416. Task (3) requires 55 selector positions with 
delayed pickup ancTsome comparing units or a five-position selector and a digit 
emitter for control. 

The problem in connection with task (1) is to get the desired zero control 
and spacing in the two cases with the same setting of the hammersplit levers (also 
called zero suppression levers) and hammerlocks, on machines where zero control 
is not performed on the control panel. On a tabulator with 89 type bars, this may 
be done as follows: 



Type bar 



Direct Wiring 



Data Cards 



Trace Cards 



Alpha. 1 -4 
5 
6-8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 



Col. 1-4, II*. 
Col. 7-9, II. 



Col. 5, II. 
Col. 6, II. 



Col. 11, I. 



Col. 12-21, II. 



Col. 22, I. 
Col. 23-29, II. 



Col. 11, I*. 

Col. 12, II. 
Emit 10. 

Col. 13, II. 

Col. 14, II. 

Col. 15, II. 
Emit 10. 

Col. 16, II. 

Col. 17, II. 

Col. 18, II. 
Emit 10. 

Col. 19, II. 

Col. 20, II. 

Col. 21, II. 
Emit 10. 

Col. 23, II. 
Emit 10. 

Col. 24, II. 

Col. 25, II. 

Col. 26, II. 
Emit 10. 

Col. 27, II. 
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Type bar 



Direct Wiring 



Data Cards 



Trace Cards 



Num. 



32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
1 
2 
3 
4 
5 
6 
7 
8-17 
18 
19 
20-29 
30 
31 
32-41 
42 
43-45 



Col. 30, II. 
Col. 31, II. 
Col. 32, II. 



Col. 41, II. 

Col. 42, II. 

Col. 43, II. 

Col. 44, I. 

Col. 45-54, II. 

Col. 55, I. 

Col. 56-65, II. 

Col. 66, I. 

Col. 67-76, II. 

Col. 77-79, II. 



Col. 33, I. 



Col. 34-40, II. 



Col. 28, II. 
Col. 29, II. 
Emit 10. 



Emit 10. 

Col. 35, II, 

Col. 36, II. 

Col. 37, II. 
Emit 10. 

Col. 38, II. 

Col. 39, II. 

Col. 40, II. 
Emit 10. 



(*) The symbol "11" denotes wiring from the second brushes (on some machines 
called "third reading") whereas "I" denotes wiring from the first brushes ("second 
reading") through a selector that separates IPs from 12<s to the X-PU of an X- 
distributor ("pilot selector") through the transfer point of which an emitted 10 goes 
to the type bar in question. 

The hammersplit (zero suppression) levers alpha. 4, 11, 25, 37 and num. 5, 
17, 29, 41 and the long hammerlocks alpha. 12, 16, 20, 24, 26, 30, 34, 38, 42 and 
num. 2 are raised. Left zero carry clips of width 3 are attached to hammersplit 
levers alpha. 6-8 and num. 43-45. 

Trace cards are distinguished by the presence of an 8 in column 80. This 
impulse is wired through a digit selector to the digit pickup of a row of selectors 
("class selectors" or "co-selectors with a controlling pilot selector") with a total 
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of 34 sets of points and to the hammerlock control hub. 

For performing task (2) in the manner specified in See. V. 2. , four external 
switches are needed. On a 416, the single-double spacing lever and the minor, 
intermediate and major control switches can be used; on machines with four 
pluggable switches there is no problem. 

Spacing every 10 lines (alternative (b), Sec. V. 2. ) may be accomplished by 
adding 1 to a counter on card cycles and using the carry (which, on a 416, is 
automatically available at the counter total exit) to initiate a minor cycle during 
which spacing takes place and the counter is cleared. 

Spacing alternatives (c) and (d) both involve inspecting a card column at two 
reading stations and taking certain action when the digit at the first brushes is less 
than the digit at the second brushes. This may be done by wiring from second 
brushes to a comparing entry and from the corresponding comparing exit to the 
immediate pickup ( rr ZFS P. U. ,c on a 416) of a selector through the transfer point 
of which the digit at first brushes is passed. In case (d), this digit is taken directly 
to cause spacing; in case (c) it is wired to the digit pickup of a selector which 
initiates spacing on the next cycle. 

Task (3) is easily accomplished if a sufficient number of selectors are avail- 
able. Since the same task can be performed on the 650 at the expense of 5 locations, 
(see Sec. XI, 1. ) the tabulator wiring will not be discussed here. 

XI. 6. SYSTEM LOADING 

The interpretive system deck, normally with 6 words to a card and an x-punch 
in 5 to get the 1 in column 6 picked up as a leading digit of the address, is loaded by 
a deck of six self-loading cards (12 in col. 1) with 7 words to a card and a card 
number in the eighth word. The first card serves the sole purpose of making a 
fixed console setting possible. The System Loading program on these cards operates 
as follows: 



8000 
1951 
0004 
0053 
0106 
0153 
0204 
0056 
0055 
0155 
0102 



70 
70 
70 
70 
70 
70 
70 
60 
11 
44 
60 



1951 
0004 
0053 
0106 
0153 
0204 
0251 
0001 
0251 
0152 
0002 




Read in the loading pro- 
gram from six load 
cards (B into 0001, C -I- 1 
into 0002). 

Read a system card (non-load), 

Go to stop if expected loc. 
B'^- loc. on card, L. 
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0007 
0205 
0057 
0051 
0052 
0103 
0156 
0206 
0054 
0203 
8002 
8003 
0151 
0201 
0207 
0101 
0105 
0107 
0104 
0157 
0005 
0003 
0152 

0001 
0002 
0006 
0202 
0154 



11 

11 

46 

21 

60 

10 

15 

10 

21 

11 

69 

24 

15 

10 

11 

44 

10 

60 

44 

60 

10 

21 

69 

00 
00 
24 
69 
00 



0001 
0252 
0152 
0080 
0006 
0251 
0202 
0252 
0070 
0252 
0)2531 

Cl] 

0154 
8001 
0070 
0105 
0070 
0080 
0157 
0001 
0252 
0001 
7777 

B 
C 1 
0000 
0253 
0001 



0000\ 
0000 ) 
0151 Y 
8003 I 

ooooy 



Go to stop if C-4- 1 <y 
B'+ n. 



Prepare accumulator for 

move. 

Store test constant. 



Move one word. 

Increase addresses by 1. 

Test for end of moving. 

Return to move another word. 
If CH-1= B<4-n, end of load- 
ing; go to console. 
Increase location by word 
count, go to read next card. 

Error stop. 



Constants 



A similar program is used for punching out the system in condensed form in 
case extensive changes, entered on self-loading, single-instruction cards, have 
been made. 

The Reset Memory Card, mentioned in Sec. V. 3. , is a load card with eight 
words. The program, essentially identical with one supplied by the IBM 650 Sales 
Research Group at Endicott, runs as follows: 



8000 


70 


1951 


1333 


1951 


69 


8000 


1953 


1953 


24 


0000 


1954 


1954 


69 


1957 


1955 


1955 


24 


0999 


1956 


1956 


61 


1958 


8003 


8003 


20 


C000Q 


0999 
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(00 



0999 

0000 70 1951 1333 



1952 8003 
0000 0000 



1952 00 0001 0000 

1957 11 1952 8003 

1958 20 0001 0999 

XI. 7. PROGRAMS 

The complete programs of the system are listed on the next 21 pages, (i) - 

(xxi). In many cases, but not always, a constant used in two programs is listed 
in both. 
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650 INTERPRETIVE SYSTEM . 

DECK CARD LOC. OR DATA INSTR. 

i « « r, ~ <- 1. GENERAL INTERPRETATION . 

1 1 1095 60 1098 1014] 

1 2 1014 10 1024 172 9V Increase i to i+1. 

1 3 1043 21 1098 8.0 lj (see cards 26 & 27 below) 

1 4 8 1 6 [ i+1 ] 110 7 Start dissection system instr. 

1 5 1107 4 6 1112 1061 Shall the loop box-be used? 

1 610613000031019 

1 7 1019 20 1023 1026 Store C, 

1 8 1026 60 8003 103 3 

1 9 1033 30 000 3 104 1 

1 10 1041 20 1045 1048 Store B, 

1 11 104 8 60 8003 1105 

1 12 1105 30 0003 1063 

4 H J £ f ~ ? 4 106 ? 1076 Is an 2 -o P 6ration called for? 

1 14 1067 10 1020 8003 

1 15 8003 6 9 8002 1081 Go to an 0| -subroutine „ 

1 17 Jni* 2L 1 1*11 10 J 6 Add the loop box to the instr. 

1 17 1046 61 8003 1061 Remove the minus sign. 

1 18 1076 60 8001 1034 

1 19 1034 10 1037 8003 

1 20 8003 65 1045 [1000*0 2 ] Go to an O^-subroutine* 

1 21 1098 60 [ i 1 110 71 

1 22 1024 01 L 

1 23 1020 69 8002 108 if 

1 24 1037 65 1045 lOOOJ 

1 2 5 1017 L^p box ( initially ) 



V Constants 



1 26 1729 69 1732 173 5") 

1 27 1 735 24 1539 104 3 

1 28 17 32 35 0002 1445 

2 1 1000 30 0003 11191 
2 2 1119 15 1030 8002 
2 3 8002 60 [ B ] 109 6 
2 4 1096 69 9 99 9 1120 



Restore the multiplication 
routine (Oj-3) to normal 
(The trig, program makes 
special use of it), 

2» TRANSFER OPERATIONS . 

0?»000 . UNC STOP 

Call outTTfor console display 

Stop. 



2 * in%o ?2 ??22 1059 -> 02-200 - CONDSTQP 

tf 6 1059 15 1162 8002L r>?rrr^!*i~n~7~Z 1 ,i -, 

2 7 8002 60 [B] 10 1 6 J console display 

2 8101601 * 1120 Stop; Go to C on PROGRAM START 

2 " 10 ^n?^ }2 }?° 9 1015 ) 2^201, TRSGN 

IS Ills Is isss "i? , £tr?f r '"™' 

2 12 1069 60 8002 1027 ft, to B if + ! 

2 13 1202 60 1009 1013 2 -202. TO EXP 
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2 
2 
2 
2 
2 

2 
2 
2 

2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 

2 
2 
2 



2 
2 
2 
2 
2 



3 
3 
3 
3 
3 



14 
1 5 
16 
17 
18 

19 
20 
21 



2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 

3 

3 
3 



37 

38 
39 
40 
41 



1 

2 
3 
4 
5 



1013 
1021 
1029 
1036 
1012 

12 3 
1027 
1035 



1 
1 
1 
1 
1 



204 
127 
040 
151 
161 
8002 
10 7 5 
1133 
1092 
1049 
1058 

1205 
1128 
8003 



3 6 1454 



114 
145 
038 
162 
030 



1 001 
1064 
1 071 
1077 
1135 



6 
30 
10 

65 
30 
69 
22 
1 5 
67 
30 
60 
15 
35 
21 

65 
35 



67 
21 
60 
60 
60 



30 
69 
30 
10 
10 



0002 
8002 
0001 
104 5 
109 5 

1023 
000 3 
103 8 



102 
000 
114 
105 
111 

[ C 

00 
800 
104 
000 

[ c 



102 3 
000 3 



00 5 
1017 
000 2 
8 001 
113 8 



1021' 

102 9' Get the exp. of prev. result 

1 3 6J 

10 12 Compare it to B; Go ahead 

112 if exp . < b, to C if exp. £ B. 



10 2 7] 2 «203 . TR 

3 5V Replace i+1 by C in 
1 4 3 J general interpretation. 



112 7^ 
104 
51 
61 



11 
11 



02-204 . TR SUBR 



800 2 
1075 
1133 
109 2 
1049 
105 8 
12 3 J 

1 1 2 8" 
800 3 
[ C ]. 

109 5 



r 



Set the O-address of the 
instruction at C equal to B; 
30 to TR to C. 



02-205 . TR OUT 

Go to C in machine language 



02*454 . N00P 

Go to general interpretation, 
107 5^| 
12 3 

1 1 7 [ Constants 
101 6 
109 6 J 

3« LOOP OPERATIONS 
10 6 4*) 02*001 * LOOP C 



107 1 
107 7 
113 5 
109 3J 



Move B to the C-address 
" position; Call out and 
increase the loop box. 



3 
3 
3 
3 
3 
3 

3 
3 
3 

3 
3 
3 

3 
3 



6 
7 
8 
9 

1 

2 
3 
4 

5 
6 
7 



1 
1 
1 
1 
1 
1 



093 
051 

54 
062 

1 55 
11 3 



18 
19 



1 100 
10 57 
1025 

1 101 
1008 
1074 

1110 
1121 



11 
46 
10 
21 
16 
2 

30 
10 

10 

30 
10 
10 

69 
30 



8 002 
105 4 
8001 
1017 
8002 
1017 

0001 
1017 
103 1 

001 
1017 
1031 

1017 
0001 



1051] Common LOOP step s 

115 5" Compare the loop box to B; 

10 6 2 if less, loop again (go to C). 

11 2 0J 

1 1 1 31 If equal, reset ihe loop box, 
1 9 5 J go to general interpretation 



10 5 7 
102 5 
109 3 

100 8 
1074 
113 5 

112 1 
102 8 



09*100 . LOOP A 



02-101 . LOOP AC 



02*110 . LOOP AB 
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3 
3 

3 
3 
3 
3 

3 
3 
3 

3 
3 

3 
3 
3 
3 
3 



20 
21 



1028 
1136 



4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 
4 
4 
4 
4 

4 
4 
4 
4 
4 
4 

4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 



22 1111 

23 117 1 
2 4 1078 
2 5 1185 

26 1010 

2 7 13 4 3 

28 1003 



29 
30 

31 
3 2 
33 
3 4 
3 5 



12 
13 
1 4 
1 5 
16 

1 7 
1 8 

1 9 

2 
21 
22 

23 
2 4 
2 5 
2 6 
27 
2 8 
2 9 
30 
31 
32 
33 



1 Oil 

118 

1138 
1031 
10 47 
1039 
1139 



10 
10 

69 
30 
10 
10 

30 
10 
10 

30 
10 



8 001 

103 9 

1017 
0001 
8001 
113 9 

000 4 

104 7 
101 7 

00 4 
113 8 



0100 

0100 
0100 



113 6 
109 3 

117 1 
107 8 

118 5 
109 3 

134 3 
1003 
109 3 

118 

134 3 

ooo n 

100 
100 
1001 j 




115 3 
1312 
1322 
1 330 
1388 



1 
1 
1 
1 
1 
1 



367 
328 
340 
363 
172 
336 



10 
30 
10 
16 
15 

10 
66 
11 
30 
11 
16 



80 2 

000 3 

JL -^ £j O 

128 3 

114 2 



1312 

132 2 

133 
138 8 
119 7^ 



02-111 . LQOP ABC 



02-010 - LOOP B 



02-011 . LOOP BC 



Constants. 



4, MOVE 
Qt-9 . MOVE 

Test N (i7e # , A);go to special 

move if N-0 

Store N-l 



Go to stop if N-l+C ;> 1000. 
Test for upward or downward move, 

r Initialize for downward move. 



102 
8 
104 
000 
132 
114 



132 
134 
136 
117 
133 
119 



8 1 







Initialize for upward move. 



1197 
1176 
8003 
8 002 
1102 



59 
66 
21 
79 
86 
52 




Store test constant * (see next p.) 
Complete initialization 

Move one word 

Test for completion 

Change addresses by 1, 
return to move next word 
If N«0, Get B and go 
to storing routine. 



~iv~ 



4 
4 
4 

4 
4 
4 



5 
5 
5 
5 
5 
5 

5 
5 
5 

5 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 

5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 



6 
6 



34 1362 
3 5 80 02 
36 1315 



37 
38 
39 



1 

2 
3 
4 
5 
6 



9 


1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 



32 
3 3 



3 
3 
3 
3 



38 
39 



1 325 
1142 
1065 



06 
07 
08 
09 
53 
00 



1009 
1023 
1045 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



50 
5 2 
55 
60 
64 
65 
67 
68 
70 
72 
74 
77 
78 
83 
285 
289 
291 
293 
294 
297 



15 106 5 
60 [ B ] 
69 222 2 

69 
24 
60 



77 
77 
77 
77 
7 7 
77 

55 
55 
55 

88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 
88 



77 
77 
77 
77 
77 
77 



O w O sJ 

5 55 5 

5 55 5 

8 8 8 8 

8 88 8 

8 8 8 8 

8 88 8 

8 88 8 

8 88 8 

8 8 8 8 

8 88 8 

8 88 8 

8 88 8 

8 8 8 8 

8 88 8 

8 8 8 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 

8 88 8 



800 2 
144 5 
109 5 



Error stop for MOVE, READ & PCH, 



!} 



80 2" 

110 2V Constants 
144 5, 



31 7 
323 
378 
381 
337 
391 
396 
66 



7 7 
77 
77 
77 
77 
77 

rj n 
11 



11 
11 
11 
11 
11 
11 
11 
11 



5. ERASABLE AND VACANT LOCATIONS , 

7 7 7 7 

111 l\ Vacant 02-codes 

7777 

777 7 

1111) 






55 5 5] 

5 5 5 5V Inter-subroutine 

555 5) 



88 8 8\ 
88 8 8 
88 8 8 
88 8 8 
88 8 8 
88 8 8 
88 8 8 
88 8 8 
88 8 8 
888 8 
888 8 
8 8 8 8 
88 8 8 
888 8 
88 8 8 
888 8 
8 8 8 8 
88 8 8 



storage 



INITIALIZATION CONSTANTS 
FOR MOVE 




Downward 


Upward 


# 


B+N-l 


-B 


## 


B 


-(B+N-l) 


4HB* 


C 


-(C+N-l) 



8 8 8 8 /-Erasable 
888 8J 



3 13 56 7 7 7777 77771 
31 1314 77 777 7 7777 



777 
77 7 
777 
777 
777 
777 
777 
777 



VVacant registers 



40 1044 77 7777 111 1J 



1 1087 69 114 119 5 



6. ADDRESS CHANGE OPERATIONS. 



10 88 69 1341 119 5 Pi*? . TR B 



-v- 



3 1089 69 1342 119 5 



6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 
6 
6 
6 



6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 
6 
6 
6 
6 



4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 



2 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 

4 

4 

4 

4 



3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 



4 5 

4 6 

4 7 

4 8 

4 9 



1195 
10 53 
1311 
1369 
8 2 
1198 
359 
117 
125 
183 
392 
199 



16 1140 

17 1341 

18 1342 

19 1072 

2 10 5 

21 10 50 

2 2 1500 



1 
1 
1 
1 
1 
1 
1 
1 
1 



6 3 
69 

7 7 

8 4 
92 



1 

2 

1 

1 

2 

3 

1 

1 

1 

8 

12 8 

108 

JL & d 

32 
144 



58 
15 
1 8 
30 
96 
74 
67 



103 
11 6 
80 
160 
104 



24 
6 5 
30 

1 5 
6 7 

3 5 



6 
3 
6 
3 
1 
4 



35 
35 
3 5 
67 



9 
2 
9 
2 
5 

9 
2 
5 





20 
60 
30 
60 
4 6 



3 5 
30 

21 



119 8 
8 02 
000 3 
1072 
[ A 1 
[1,4,7] 
8002 
000 3 
8 00 
00 
104 
120 




0l»8 . TR C 

Common stegs. 

Set amount of shift for TR A, BorC. 

Get the instruction ("K) 
located at A. 



Separate out its A-, B-, 
or C-address. 



01 
000 4 
00 7 



Compare this to B. 
Go on if equal, to C 



Constants 



if unequal. 



65 8 00 3 
65 110 3 



135 91 
135 9 
135 9 
119 8, 



116 3 02-005 . SET C 



116 3 



65 11041163 



1116 

127 4 
108 

128 9 
104 
00 
116 
12 6 
107 

[0,3,6] 

126 5 
8 
00 
8 00 
111 
10 2 
102 

12 6 
[0,3,6] 
[ B ] 



5 
3 

7 
3 

] 



3 
3 
3 
5 
3 
3 
3 
5 



1 
1 
1 
1 
1 
1 
1 
1 



2 

1 

1 

2 

3 

1 

1 

1 

800 
12 8 

1 

JL fo 



1 
1 
1 
1 
1 
1 

2 
2 




} 



02-050 . SET B 
Op-500 . SET A 
Common steps 
Set amounts of shift. 



\ Set the address for storing 

J the modified instr. at B # 

Get the instruction (1) 

located at B. 



9"\ 
7 

4 

2 

91 

7 

4 



2 

9. 

8' 

2} Store the right end of B. 

2 

4> Destroy the old A~, B- or 

8 J C-address 

8] 

\ Replace it by C 



si 

i) 

iy Attach the right end; 
5J store at B. 



3 



6 



12 6 7 
110 8 
109 5 



Constants 



50 1006 651023 1079 0g«006 . ADD C 

51 1007 6 6 1023 1079 0g»007 . SUB C 



-vi« 



6 
6 
6 
6 
6 
6 

6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 

6 

6 
6 



52 1079 30 0007 1004 Common to both 

53 1060 65 1023 1129 2 -060 . ADD B 

54 1070 6 6 10 2 3 1129 02«070 . SUB B 

55 1129 30 0004 1004 Common to both 

56 1600 65 1023 1189 02*600 . ADD A 

57 17 00 6 6 1023 1189 02»700 . SUB A 

58 1189 30 0001 1004 Common to both 



68 

6 9 

70 
71 

7 2 



1004 
1068 
1399 
13 57 
1313 
1370 
8 2 
1348 
1106 
1156 
1267 

1360 
1123 
-107 3 



20 
65 
30 

69 

JL \J 

60 
4 6 
10 
1 1 
21 

21 
60 
60 



1260 
104 5 
00 
136 
126 
112 



115 6 



15 
26 
26 

B ] 



106 8" 
139 9 
135 7 
1313 
137 
800 2 
13 4 8 
110 6 
12 6 7 
12 6 7 
109 5 

109 5' 
134 8 
128 9 



Common to all A0p & SUB o#, 
Store shifted ±c"; ~~ 
Set address for storing 
modified instr* at g 

Get I 



Add ±C to J B|j store 
(with original sign) at B, 



Constants 



7 
7 
7 
7 
7 
7 



7 
7 
7 
7 
7 
7 
7 



1 8000 70 1951 



2 
3 
4 
5 
6 
7 




1 
2 

3 
4 



33 
24 
94 
78 
75 
44 



60 
24 
20 
24 
2 4 
60 



15 1038 60 



2 4 2 
27 8 

017 
98 6 
195 2 



8 1042 45 1364 



9 1377 69195 



198 
197 
106 
19 5 
10 3 



133 
112 
119 
117 

134 
10 4 
1 37 



7. PROGRAM LOADING , 

Read the first program card. 

Set an instruction for return 

from the READ routine 

Reset 000 

Reset the Loop box 

Reset the card counter 

Test the word count, n. If / 0, 

go to READ to move n words 



17 1242 70 1951 1344 



1 3 7 6\ If n»*0, put the problem no, 

1 3 9 7 J into the punch band, 

13 4 5") decide if tracing is required 

1 1 8 J (dep, on contents of 1976*), 

11261 store the address of the 

17 2 9 V first programmer instruction, 

go to execute it (gen, int.) 
1 1 7 J |3o 0003 1019 if not tracing 
121 1980 1386 if tracing 

Constant (return from READ) 



8 
8 


1 
2 


1401 
1131 


01 
60 


8 


113 1 

14 4 5 


8 
8 


3 
4 


14 
1174 


69 

24 


1318 
12 7 8 


117 4 
14 8 6 



8. READING OPERATIONS 

1 Og'AOl . READ CONS 

i Stop* On start P read console 

switches, go to storing routine* 
1 02=400 . READ 
r Set return instructions, go to 



-vii- 



8 

8 
8 
8 
8 
8 
8 
8 
8 
8 
8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 
8 
8 
8 



5 1486 69 1241 116 8j Steps common with PCH 

Return from common steps. 
Read a card. 

Go to stop if loc. L on the card 
differs from progr. first loc* 
(the stop is in MOVE) 
Go to stop if B»+n > C+l 
(n=word count, B 1 "current 
first location) 



Calculate and store 
test constant 

Complete initialization 

Move one word from the read 
band into place 
Increase addresses by 1. 

Test for end of moving « 

Return to move next word 

If B«+n » C+l, end of READ 

Add n to B* to get new B 1 , 
go to read next card« 

Constants 




33 1188 69 

3 4 1318 60 

3 5 1241 21 

36 1324 24 



195 3 

128 5 
127 7 



80 

173 6 
132 



3 1 
V 



1 1410 6 



196 4 116 8 



9 
9 
9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



2 
3 
4 
5 
6 
7 



168 
81 
2 32 
36 1 
319 
134 



24 
65 
10 
30 

1 5 

2 



1 <i f <J 

10 2 3 

104 5 

000 3 

10 2 4 

129 1 



9. PUNCHING . 
02«U0 . PCH 




2 



108 11 Common w ith READ 

12 3 2 

1 3 6 1 1 Prepare for testing 

1 3 1 9 f and initialization. 
113 4 

JL <z> f & , 



PCH program only 

Go to stop (in MOVE) if 

B > C+l. 

Set the word count 

Let B 1 -first loc. not yet 

punched. If B»+5 < C+l, 

go to punch 5 to a"~card. 

If not, go to further testing. 

Prepare to move n words 
to the punch band 



-viii- 

9 2 1 1950 15 197 7 1987 

9 2 2 1987 10 1990 8002 

9 2 3 8002 6 9 [B»] 8003V Move one word 

9 2 4 8003 2 4 [19 7 9] [1990+n] j 

9 25 1991 60 1986 1971 Go to punch 

9 26 1992 15 1024 198 7) 

9 27 1993 15 1024 1987 

9 28 1994 15 1024 1987}- Return to move 

9 29 1995 15 1024 1987 another word 

9 30 1996 15 1024 198 7J 

9 3 1 19 7 1 10 10 2 4 1 9 6 7l Increase the card no. 

9 32 1967 21 1986 1 9 9 7J 

9 33 1997 71 19771947 Punch a card 

S 34 1947 60 1977 194 6) 

S 3 5 19 4 6 10 19 7 8 1 9 6 5 V Set the location no. 

S 3 6 19 6 5 2 1 19 7 7 1 9 6 6 J for the next card 

9 37 1975 10 8001 19 4 91 If B»-C+l, punching 

9 38 1949 4 4 19 7 3 1 9 5J is completed. 

[This const* keeps trace orders 

9 39 1976 30 0003 101 91 {inoperative during non-tracing. 

9 40 1964 21 1977 198 9 (This const, will be 00 0000 1095 

9 411969 0005 {.during tracing (PCH inoperative) 

9 421990 99 9 9> Constants 

9 43 1963 24 1978 199 1/ 

9 44 1974 69 800 3 (When the PCH progr. is loaded 

9 45 1958 24 1978 199 ljjwith program loading, this instr. 

twill load into 1963 (I). 
10. TRACING 

10 1 1061 21 1980 1386 Expansion of gen. int. 

10 2 1386 30 0003 1946 Store modified instr. for trac. 

10 3 1946 201023 1947 Store C 
10 
10 

10 6 1949 20 1045 1950 Store B 
10 
10 
10 
10 

10 11 1965 6 5 8 002 1 9 6 6\- Call out A if Oi / 
10 
10 
10 

1Q 15 1968 21 1 9 8 £ 1969 Store A or for tracing 




10 

10 

10 18 1971 15 1999 8002V Store B for tracing 

10 

10 

10 2 1 1973 6 9 19741 Store prev. result for tracing 

10 2 2 

10 23 1975 69 1017 198 7} Store loop box for tracing 

10 2 4 

10 25 1988 65 1098 198 91 Store instr. no. for tracing 

10 2 6 

10 2 7 

10 28 8002 6 9 f i+1 ] 1 9 9 1 V store original instr. for tracing 

10 29 1991 

10 30 1992 






10 
10 
10 
10 

10 
10 

10 
10 

10 
10 
10 
10 
10 
10 

10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 



11 

11 
1 1 

1 1 

11 
11 

11 



1 2 
12 
12 
12 



3 1 

3 2 

3 3 

3 4 

3 5 
3 6 

37 

3 8 

3 9 
40 
41 

4 2 
4 3 
4 4 



4 
4 
4 
4 
4 
5 
5 
5 
5 
5 
5 



1 

2 
3 
4 
5 
6 
7 
8 



1 

2 
3 
4 



6 
7 
8 
9 




1993 
1 994 
1995 
199 6 



4 50 

331 

4 51 
3 5 8 



14 52 
1332 
1091 
1165 
1334 
8 2 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



14 4 5 
1404 

1414 
1427 
14 30 
14 3 9 
14 4 7 



8 02 



82 
49 2 
4 49 
4 57 



15 
20 
71 
60 

69 
24 

69 
2 4 

69 
24 
6 5 
16 
69 
2 4 

09 
6 
6 9 
04 
36 
30 



1 
9 
9 
9 

9 58 4 4 



21 



9 1443 24 



2 4 
65 
30 
1 5 



102 4 
198 6 
197 7 
126 

197 6 
1061 

115 4 
106 1 

197 6 
1061 
109 8 
117 3 
1141 

[i + l J 



5 400 

00 3 
6 

198 
01 
196 5 



100 9 

102 3 

14 3 

3 

144 3 



4 [ C ] 



12 8 9 
104 5 
3 
146 



199 4V 
199 5) 
199 6 
110 5 



133 
109 

135 
109 



I) 



Store card no. for trace cards 

Punch a trace card 

Return to normal gen, int. 

02*450 . START TR 

Modify gen. int. to include 

tracing expansion 

02"451 . STOP TR 

Restore gen. int. to normal 



10 8 3 
154 2 
1599 
120 7 
8 02 



11 14 7 3 

12 1483 

13 14 93 



2 4 
65 
30 
1 5 
61 

21 
6 5 
4 5' 



128 9 
104 5 
3 
1310 

[ B ] 

127 8 

128 9 
149 6 



13 3 2] 02-452 . ST TR ERAS 

10 9 1} Modify gen, int. to include 

116 5) tracing expansion. 

1 3 3 4 ( 

8 2 V Replace the present progr. 

1 9 5 J instr. by a NOOP 



Constants 



Word count ^ trace identif . 

Make PCH inoperative when tracing. 

Will be 30 0003 1019 when not taac. 
'When the trace program is 
19 6 8j-| loaded with program loading, 
.this will load into 1963 (l). 

11. STORING THE RESULT 

(Common to all math, routines ) 



\ 



Store C in 000 for prev. res. and 
in 1009 for cond c transfer ops. 
If C-000, go directly to gen. int. 

If C ^ 000, store C irt C, 
go to gen. int. 



10 9 5 Constant 



14 9 2 
144 9' 
14 5 7 
800 2 



12. ADDITION AND SUBTRACTION 
pl-1 . ADD 



Get B 



5 8002 60 [ B] 147 3. 



154 2 
15 9 9" 
120 7 
800 2 
14 7 3, 

148 3 

149 3 
149 7 



SUB 



Get -B 



Common 



If A*000, get -K directly. 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 



1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 



40 
41 
4 2 
4 3 
44 
45 
4 6 
47 
4 8 
4 9 
50 
1 
2 
3 
4 
5 
6 
7 
8 



5 
5 
5 
5 
5 
5 
5 
5 

5 9 

6 



12 
1 2 

12 



6 6 
6 7 
6 8 
6 9 
70 



149 
14 
8 
14 4 

4 

4 

4 

4 

4 

5 

5 

6 

6 

6 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



30 
1 5 



89 
9 6 



8 
29 



4 8 

5 4 



6 

5 

4 

4 

4 

4 

4 69 

472 

41 

3 

1 

2 



1 

3 



2 
9 
3 
3 
1 
9 
2 



13 10 

14 6 
14 9 
1433 
1456 



00 3 
140 9 
[ A ] 



60 
30 
20 
21 



60 
36 



1 
6 
3 
4 

6 
3 
4 
4 
6 
6 
4 
6 
6 
3 
1 
6 
2 
6 
6 



61 
6 
60 
31 
01 



1 
1 
1 
1 



00 
126 

1 27 
1 27 

00 

1 29 
1 25 
8 00 
126 
00 4 
143 3 

63 9 
69 7 
25 
272 
2 

148 9 
00 3 
1652 
14 56 
12 6 5 
2 
8 02 

145 6 
126 5 

149 4 
8 00 2 

8 02 
800 3 
002 

146 1 



000 
141 
142 
8 00 
127 
147 
125 
127 

C 

8 

129 
126 
127 
125 



5 
8 
2 
1 
8 
2 

2 

] 
1 

3 

5 

2 





Get A if A / 000 



A - Ai ,a # Split up 
and store 

B - Bi ,b. Split up 
and store 



Put |b-a| into shift instr. 



-If a*b, form Ai+Bi*c] 



0|lf |c| | > 10, C, « c{Ao. 

If Ci < 0, get c*a+l f 
combine with Ci , 
go to store. 

Same if C] > 

If Ci^O, go to store machine 0, 



If jd| < 10, get -C, by 
shirting; correct exp. 

If a / b, prepare to shift « 




Ifjb-a| > 10, get OA or 
OB ? go to store. 



If b < a, shift Bt . 



If b > a, interchange 

Ai and B^ replace a by b, 



Constants 



13 
13 
13 
13 
13 

13 
13 
13 
13 
13 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



3 
3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 



1 
2 
3 
4 
5 

6 
7 
8 
9 
10 



1 

1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2 
2 
2 
2 
2 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 



2 8 
2 9 



3 
3 
3 
3 
3 
3 

3 
3 
3 
3 
4 
4 
4 
4 
4 
4 
4 
4 
4 8 

4 9 
50 

5 1 
5 2 
5 3 




1 
2 
3 
4 
5 

6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 



10 84 
159 2 
15 49 

15 57 
8 02 

10 86 
1692 
1649 

16 57 
8 02 



5 78 

5 86 

59 

59 

60 
8 
159 



1 

1 
1 
1 

1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



6 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 
5 



5 
1 

1 

2 

3 

3 

4 

4 

5 

6 

7 

3 

5 88 
59 1 
54 6 

5 6 4 

6 4 8 
2 59 



3 
6 

5 
2 
7 
5 
1 
5 
3 

7 
3 
7 
6 
5 
3 
1 



717 
273 
677 
533 
628 

4 8 4 
7 43 

9 9 
54 
7 
5 8 
53 
63 
571 

5 28 
5 81 
5 32 
5 39 



24 
65 
30 
1 5 
60 

24 
65 
30 
1 5 
61 



2 

6 

4 

3 

1 

6 

6 

3 

4 

2 

2 

6 

3 

4 

2 

6 

1 

1 

2 



1 

5 

5 



5 







4 

1 







4 

1 

8 

5 

8 





4 6 
4 4 
60 
19 
3 5 
6 5 



3 

4 
6 
3 
6 
3 
3 
4 
.4 
1 
1 



5 
4 
9 
1 
9 
1 
5 
4 
6 
5 
6 

3 5 
60 
30 

4 6 
18 
1 7 
3 5 



128 9 
104 5 
00 3 
156 
[ B ] 

128 9 
104 5 
3 
166 
[ B ] 

128 3 
128 9 
159 6 
3 
150 8 

[ A 1 

002 
1515 
1270 

127 7 

128 3 
2 
154 7 

JL £$ O O 

8 2 
156 8 

27 

2 8 

59 

59 

2 5 



1 
1 
1 
1 
1 

127 
4 
8 3 

01 
167 7 
15 8 
00 3 
128 
000 


125 

1 50 
8 
8 
00 
8 2 
2 
158 1 

128 5 
128 5 
00 2 



1 
2 
2 
3 
7 
1 
1 
7 
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13* MULTIPLICATION 

15 9 2 0^3 . MPY 
1549] -" , 

1 5 5 7 V Get B 
800 2 
1 5 7 8 J 

1 69 2 0i"5 . NGMPY 

16 4 9] 

1 6 5 7 V Get -B 
800 2 
15 7 8 J 

15 8 6 Common steps 
-Get A 

If AtQ, CO directly* 
A*Ai,a. Split up and store 
B*Bi,b. Split up and store, 



1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



Calculate and store 
49 - (a+b) 

Go to error stop if exp* of 
prod, would be out of range 

Ci - At • B) 



Common with DIV . 
If C{ > 10, C, - CJ/10. 
41 If Cj < 10, prepare to 

91 

4 

8 

3 

3 

8 

1 

8 

2 

9 

9 

5 



correct exponent 

| If Ci«10 (due to carry in 
.rounding), go to correct • 



r Calculate c « 
attach it to 



s a+b-49-corr. j 



J 



-xLi- 



13 54 1595 6 9 5555 1338 Error stop, for MPY and DIV. 



14 

14 

14 

14 

14 

14 

14 

14 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



14 

14 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 

1 
1 
1 



4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 
4 

4 
4 

4 

4 
4 
4 



15 
15 
15 
15 



61 
62 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

22 

2 3 



2 

2 
2 
2 
2 
2 
3 
3 
3 



3 3 
34 
3 5 



1 
2 
3 
4 



660 
560 
508 
568 
580 
281 



1338 
1316 

10 85 
1642 
1699 
1607 
80 02 
1626 
636 



1 

1 

1 

1 

8 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



43 
46 
56 
02 
47 
606 
613 
6 20 
680 
637 
693 
650 
407 
6 
7 
2 
8 
4 
5 
1 



6 2 
62 
67 
63 
64 



1610 
1659 
1654 



3 6 167 3 
3 7 15 2 7 
3 8 1374 



1350 
1169 
8 02 
1137 



61 
60 
60 
49 



65 

60 

2 4 

65 

30 

1 5 

60 

21 

65 

45 

30 

15 

60 

60 

30 

20 

31 

60 

30 

20 

60 

35 

2 

6 

3 

2 

6 

1 

1 

4 

4 



1 

5 
1 
7 
6 
8 
6 
4 
20 
60 
64 

60 
60 
50 

4 5 
45 
69 



30 
15 
67 
60 



8 2 

80 03 



1 2 8 
104 
000 
1 61 
[ B 
128 
28 
64 
00 
65 
A 



1 
1 

1 

r 



00 2 
1267 

1 27 
128 
000 
1 29 
8 00 

000 2 
126 

1 27 

129 

29 
6 5 
26 
62 
59 
28 
29 
26 



1 
1 
1 
1 
1 
1 
1 
1 



162 6 
168 
333 3 



3 

10 2 3 

[ B ] 

8 00 2 



15 7 8") 
157 8 

16 5 5V Constants 



000 lJ 



144 

15 3 



Go to store after error stop 

1U. DIVISION . 
0i-4 . DIV 

Get B 

BO tested below 




0h 



B«Bi »b» Split up and store 



Store shifted Aj 



Calculate and store b-a~50. 



\Go to stop (in MPS) if exp. of 
J quotient would be out of range, 

ol Ci«Ai/Bi; go to exp. corr. 

7 J (in MP?) and storing routines* 

16 7 3] 

16 6V Constants 



137 4 
144 5 
133 8 



1169 
800 2 
113 7 

144 5 



Go to stop if B*0, 

Set C-0 directly if A^O 

Error stop for B*0. 

l fi. ABSOLUTE VALUE . 

Par22g> M 

Get |b| * C, go to store it. 



-jdii- 



15 



5 10 22 6 7 



113 7 Constant 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6. 

6 

6 

6 

6 

6 

6 



16 
16 
16 
1 6 
16 
16 



1 

2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 



2 7 

2 8 

2 9 

3 
3 1 



1300 
1214 
8 
120 
67 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 
1 
1 
1 
1 



48 
29 
51 
35 
45 
52 
9 



20 
21 
2 2 
23 
21 
22 
1 8 
23 

1 9 

2 4 
210 
2 28 
2 44 
2 37 

IZi & & 

20 8 
215 

2 

3 

9 

4 

1 

2 

4 

1 

2 

5 

2 



4 
3 
2 
6 
2 

3 
3 
1 
7 
6 



21 
2 3 
24 

5 
41 

1 5 



30 
15 
60 
46 
69 
60 
30 



4 
2 
6 
3 

2 
6 
1 
1 
3 
1 
1 
2 
6 
2 
1 
3 
4 
6 
1 
1 
1 
1 
3 
6 
1 
3 
1 
2 
6 
4 



60 
19 



60 



1670 



000 
121 
B 
67 
4 44 
800 
000 
122 
127 
8 00 
000 
125 
123 
125 
1 24 
00 
125 
8 
127 
8 
129 
127 
00 
123 
8 
129 
10 5 
8 2 
8 001 

105 
127 

141 
127 
8 
124 
115 
127 



8 00 2 
2 
127 8 
1270 
2 



16. SQUARE ROOT 
02-300. SORT 

Get B 

Stop if B < 0. 

If B«0, C«0 directly, 
B*Bi,b. Split up and store 

x *l+ B 22Bi 



r Calc. and store Bj/xjj. 
If (Bt/j^ - xj < 10~ 8 , go to end. 



x^ « 1/2 (B]/^ + x n ) 



* Calc. 1/2 (b+50). 






If b is even, go to end. 

If b ia odd, Ci - x» vlo. 

Equip with exp. , go to store 
If b is even, Cj-x N . 

y Constants 



03 1622 



17 
17 



1 1301 69 



1514 1467 



17* EXPONENTIAL 
02*301 . ME E 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 



2 1351 69 1464 1467 Q P -351 . EXP 10 



3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 
6 
7 
8 
9 



20 
21 



2 8 
2 9 



4 2 
43 



4 8 
4 9 



14 67 
1425 

14 88 
80 02 

15 53 



5 6 




Common steps . 

Set an instr. for EXP E or EXP 10. 

Get B*Bi ,b. Split up 
and store. 



Set an instr. distinguishing 
between pos. and neg. exp* 



Form b-49. Go to special 
routines if |b-49| 2 10 * 



Set instr. for left or right 
shift of |b-49|. 
Go to EXP E routine 
Go to EXP 10 routine 

For eB-lO 1110 , calc. |mBi [• 

For 10B, get |B] |. 

Shift |B|»n+\. (0 £ \ < 1). 
Store \. 
Store n. 

If |I| > 50, go to stop* 



10^ * [l+a^X+.^tayX 7 ] 2 



-XV- 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

17 
17 
17 
17 
17 
17 
17 



7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 



1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 

7 



17 



18 
18 



5 

5 

5 

6 

6 

6 

6 

6 

6 

6 

6 

6 

6 

7 

7 

7 

7 

7 

7 

7 



7 

8 

9 



1 

2 

3 

4 

5 

6 

7 

8 

9 



1 

3 

3 

4 

5 

6 



7 7 
7 8 

7 9 
80 
81 

8 2 
83 
8 4 
8 5 
86 
8 7 
8 8 
8 9 



3 
4 
5 
6 
7 
8 
9 

1 
2 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



480 
687 
498 
66 
25 
8 3 
80 
38 
22 
679 
470 
428 
632 
641 
604 
616 
675 
431 
462 
615 

514 

464 

4 41 

52 

42 

4 4 

8 

8 



8 

4 8 
6 
58 
12 
6 



4 

4 

5 

5 

5 

5 

6 

6 

5 

51 

46 

66 

43 



4 
4 
4 
1 
2 
4 
7 



9 6 1393 



1 1302 



65 
31 

3 5 
60 
10 
10 
21 
60 
30 
64 
35 

4 4 
1 
5 
5 

1 

6 



60 

22 
22 
60 
65 
21 
4 9 
30 
35 
43 



02 
06 
11 
10 
10 



8 
000 


8 
127 

1 58 
128 
166 
000 2 
128 5 
01 
1431 
0002 
144 4 
00 02 
8 02 

127 8 
1434 
1615 
143 4 

149 5 
149 5 

8 00 3 

128 5 



429 4 



09 
025 

1 74 
7 29 
543 
627 
5 12 



For B > 0, round l6\ 
equip with exponent n+50, 
go to store. 



r\ 



r For B < 0, calculate 10 



Equip 10 ^ with exponent 

49-n (if I0"k < 1) or 

50 (if 10-A - 1); go to store. 



Go to stop if b > 59, 
set C«l if b £ 39. 



15 5 1 
150 2 
1553 
168 7 
143 8 



146 
146 
481 
005 
264 
491 
X. JL f~> 

173 
35 7 

8 8 4 
277 6 






Constants 



5 



oj 



69 6 66 6 1338 Error stop for exp. , log, & trig. 

lBo LOGARITHM , 
69 1807 1818 Og-302. LOG E 



2 1352 6 9 1815 1818 0g«352 . LOG 10 

Conation steps » 



181 
18 2 
183 
8 



8 



24 
30 
15 
60 
30 
20 
4 4 
4 6 
60 
11 



1821 
00 3 

183 6 
C B ] 
2 
129 1 

184 7 
139 3 
8 3 
1812 




Get B-Bi ,b # Split up, 
store exponent 

Go to stop (in EXP) if Bj £ 0. 



-xvi- 



1 8 

18 

18 

1 8 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 

8 
8 
8 
8 
8 
8 
8 
8 
8 
8 

8 
8 
8 
8 
8 
8 
8 
8 
8 
8 




1 8 7 1 V If B, < e, set u * B,//e 



v«,5. 



If Bi > e, 



set u • Bi/e 1 * 65 



u+1 



log- t~ + * 



t[2 



1-t 
f 



¥ 



2 + io 



n 



t 10 ]+v. 



End of common steps 



For LOG E, store L, 
calculate Mb-, . (bf 



■b-50) 



8 82 
8 40 
8 10 
819 
8 80 



35 

46 
49 
6 2 
68 



138 
000 
8 
00 
183 
00 
184 



86 
27 



8 9 5 
8 2 2 
8 8 2J 

8 4 0^ 
810 
819 
8 8 
89 3 
8 4 
44 
8 6 
8 6 
87 



For LOG 10, store mL, 
calculate bj* 



Common steps 



•(Go to stop on loss of more than 
(.two digits of accuracy 
JGo to store machine zero 
(if C has 3even zeros. 



-XV11- 



18 
18 

18 
1 8 
18 
18 
1 8 
18 
18 



8 



1 

1 8 
1 8 
18 



18 

18 
8 
8 
8 
8 
8 
8 
8 



1 

1 

1 

1 

1 

1 

1 

18 

18 

18 

18 

18 

18 

18 

18 

18 

18 

18 

18 

18 

18 



19 
19 
19 

19 
19 
19 

19 
19 
19 



6 9 
70 

71 
72 

7 3 
74 
7 5 
7 6 
7 7 

7 8 
7 9 
81 
82 



18 75 
1834 
18 43 
804 
811 
8 78 
8 08 
8 79 
8 87 



65 
31 



3 
4 
1 
1 
1 
1 
6 



5 
6 
6 
5 
5 
6 




17 50 

18 54 
1874 
18 72 



4 4 
30 
1 5 
1 6 



83 1893 01 



84 
8 5 
86 
8 7 
8 8 

8 9 
80 

9 



9 

9 

9 

9 

9 

9 

9 

9 8 

9 9 

00 

01 

2 

3 

4 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



8 



8 
8 
8 
8 
8 
8 
8 
8 



7 
5 
6 
2 
8 
7 
3 
6 

6 
9 



1 

3 

1 

2 

6 

6 

1 

3 
8 8 
86 
832 
8 48 
806 
8 42 
776 
8 56 
8 55 
8 53 
8 52 
8 96 
8 6 5 



1 1303 

2 1614 

3 1721 



4 

5 
6 

7 
8 
9 
10 
11 
1 2 
1 3 
1 4 
1 5 



1304 
1664 
1622 

15 79 
15 36 
15 45 

14 55 
80 02 
157 4 

15 8 2 

16 89 
12 67 



10 
10 
60 

60 

5 
46 

1 6 
1 9 
10 
20 

50 
23 
43 
18 
22 
28 
40 
66 



69 
24 
69 

69 
24 
69 

24 
65 
30 
15 
60 
21 
30 
2 1 

C ] 



8 00 3 
2 
002 
180 8 
1272 
1832 
1272 
1832 
8 2 

185 4 
001 
189 6 
189 6 



1 
1 

2 
6 



27 
27 



71 8 
530 



1872 
5000 
204 9 



25 8 
4 29 4 
1818 
2 222 
5714 

6 66 6 



125 6 

126 7 
152 6 

170 7 

126 7 
157 6 

128 3 
8 2 
00 3 
156 1 

C B 2 

127 8 

0002 
1252 
[ ] 




L+Mb T for LOG E 
.mL+b, for LOG 10 



Normalize and round* 



180 41 
1863 
180 4 
180 4 

1835 



182 
188 
183 
281 
659 



9) 

1 

1 

8 

7 



r In case a carry occurs 
in rounding, correct,* 



Stop on loss of two digits, 
return to log, progr. on start. 



1874 
908 6 

0052 

509 3 
4819 
1818 
222 2 
285 7 

666 6 
000 1 



1614 
172 1 
157 9 

166 4 
162 2 

157 9 

153 6 

154 5 

14 5 51 
80 2 

15 7 4 

158 2 

16 8 9 
12 6 7 

C 3 



*Card out of place. 



> 



Constants 



19, , SINE AND COSINE , 
02*303 . SIN R 



22224* COS R 



Common steps . 

(this step needed in degree prog.) 



Get B*B T ,b. Split 
up and store. 



xviii- 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 
9 
9 
9 

9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 



19 
19 



1 6 
17 
1 8 

1 9 
3 
31 

2 2 
2 3 



2 
2 
2 
2 
2 
2 
3 
3 
3 
3 
3 
3 
3 
3 
3 
3 



40 
41 



2 
3 
4 
5 
6 

8 
9 

1 
3 
3 
4 
5 
6 
7 
8 
9 

1 
2 
3 
4 
5 



6 6 
6 7 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



75 
38 
95 
98 
49 
08 
66 
74 



60 
68 



27 

82 
88 
82 
38 



28 

7 

7 



42 
47 





6 

6 

7 

7 



4 

4 

3 

20 

24 

88 

44 

45 
712 
668 

3 

8 

9 

6 

2 

9 

7 



1 
3 
4 
3 
4 
6 
2 
6 
3 
1 
4 

6 
1 
3 



46 
10 
61 
10 
10 



20 
60 



80 



8 
169 
29 
56 
25 
71 
56 
27 
27 
44 
00 
127 

6 
6 
2 
5 
2 



1 
1 
1 
1 
1 

1 
1 




1 
1 
1 
1 
1 
1 
1 



70 
37 



2 

1 

8 

9 

2 

1 

9 

8 

2 

4 

1 

6 

5 

3 

5 

3 

5 

3 

2 

5 



1618 
1252 
2 



y »i ] 



147 
188 
8 00 
800 
8 00 
188 
8 00 
170 
8 
126 
8 00 
1 67 
8 00 
8 00 
800 
000 
800 
138 
128 
8 00 
8 00 
125 
8 3 
157 2 
8 00 

1 49 
125 
8 00 
1 53 
125 
8 00 
1 68 



6 
6 
3 
3 
3 
6 
1 
5 
1 
6 
1 
6 
1 
1 
3 
9 
1 
3 
5 
1 
1 



53 
29 
7 4 



7 

6 
7 
4 
2 




6 

7 

6 

6 

4 

6 

6 

70 
171 
172 
128 
168 
171 
17 3 
174 

7 

7 

7 

4 

7 

2 

7 

7 



1 

1 
1 
1 
1 
1 
1 
1 



17 
17 



28 
7 4 



4 

6 
4 
7 
3 
4 
4 
4 
2 
4 
2 



7 

7 

2 

7 

2 

6 

5 

5 

7 

7 

7 

7 

28 

74 

74 
171 
166 

7 

3 

2 

3 

7 

2 

3 

2 

3 



8 A Tests in sine routine: 

5 

9 V Test 47-b. 

8 1 

6>If b-47, test |B|-.0025 

5 If |b| < .0025, go to store C-B 

1 r If |B| > O 0025, test 49-b, 
9 " go to special routines if 

4 J |49-.b| > 10. 

7") ' ■ - 

2 

8 1 For |49-bj < 10, set 

5( shift instruction, 

8 
6J 

0"| If b > 52, stop to indicate 

5 V loss of accuracy; continue 
5 J on start 

61 

8 

3 

6 

3 

1 

1 

4J 

3la - 4G for SIN, 1+4G for COS. 

2J 

71 If G £ 1, a ! *a 

21 if l< a £ 2, a'-2-a 

8J 

lllf 2 < a£ 3, a«-2-a 

If 3 < a< 5, a f -a-4. 

°1 

.These angle reductions 

carried out. 



Compute the fractional 
y part F of I/2nr, then 
G»F if F ^ 0, G-F+l if F < 0. 



1 
1 
1 
1 
1 
1 
1 
1 
1 



8 



8J 

4llf a'-l, s«t S-l. 

8J 

41 

5 

2 

8 



7 

6 

8 

8 



1 



9 



^S"'C 1 a l +C3a , 3+ . •*+C 9 a* 9 . 



-xix- 



19 
19 
19 
19 
19 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

19 

19 

19 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



19 
19 
19 
19 
19 

19 
19 
19 
19 



20 



74 
7 5 
7 6 
7 7 
7 8 
7 9 
80 
81 



8 8 

8 9 

9 
91 
9 2 
9 3 

9 4 
9 5 
9 6 
9 7 
9 8 
9 9 
00 
01 
02 
03 
4 
05 
06 
07 
8 
9 
10 
11 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



3 89 
275 
7 34 
94 
26 
4 
72 

7 9 

8 3 
7 22 
284 
395 

79 
37 
63 
25 
87 
86 
76 
80 



2 56 

526 



7 

6 

9 

3 

3 



7 
6 1 



7 
6 
1 
1 
1 
5 
5 
2 
8 
2 

4 
4 
8 

4 



19 

30 

60 

10 

19 

65 

35 

4 4 

31 

35 

31 

3 5 

60 

36 

11 

46 

10 

11 

46 

60 

20 
69 
68 
11 
60 
47 
99 
30 
35 

15 



07 

64 
15 

49 



125 
00 
8 
139 
128 
8 00 
00 
138 
000 
002 
002 
001 
8 002 



8 002 
128 6 
1690 
169 
139 3 
186 5 

1272 
188 6 
8 2 
188 6 



7 500 



00 
915 

1 

6 

6 

9 





1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 



2 
7 
3 
7 
6 
6 
2 
2 
7 
3 
3 
3 
7 



5 

4 

4 

6 



2J 

9^ 



7 

3 

9 

2 

4 

7 

7 

8 

2 

7 

9 

7 

3 
2 6 
72 
28 
4 4*5 
4 4 5) 
380 
28 3 



r 



Round and normalize, 
go to store 



Go to stop 
set a*0 if 



if b ^ 59, 
b£39. 



157 5 
174 2 
162 7 
174 2 
1574 



1716 

17 16V Constants 



943 
484 
765 
679 
711 
631 
005 



19 112 1353 69 1384 1739 
19 113 1354 69 1385 1739 



114 
115 
116 
117 
118 

119 
120 
121 
12 2 



1 
1 
1 
1 
1 



739 
50 5 
585 
698 
5 55 



13 8 4 
138 5 
126 1 
1167 



24 
69 
2 4 
20 
6 

21 
21 
35 
1 7 



129 3 
1261 
153 9 
128 9 
116 7 



002 

453 2 



150 5) 
158 5 
1 69 8J 
155 5) 
157 8j 

130 31 
130 4 
129 3 
9 3 4 8 



02»353 . SIN D 

Og-354 . COS D 

Common steps 

Set an instr. at the end of 
MFY to get out to SIN R or COS 
Go to MP5T with & in place of A 
and 2rr/360 in place of B. 



r Constants 



R. 



1 1305 69175 8 1761 



20. ARC TANGENT . 
P2 "?0? » ARC TAN R 



~XX« 



20 



2 1355 69 1908 1761 Op«355 . ARC TAN D 



Common steps 



20 


3 


1761 


24 


126 4 


1917 


20 


4 


1917 


30 


00 3 


192 5] 


20 


5 


19 25 


1 5 


192 8 


800 2 


20 


6 


80 02 


60 


129 7 


194 1 


20 


7 


1941 


21 


190 


20 


8 


1900 


3 


00 2 


175 7 


20 


9 


17 57 


21 


126 5 


176 8. 


20 


10 


1768 


46 


1771 


17 7 21 


20 


1 1 


17 71 


69 


177 4 


177 8< 


20 


1 2 


1772 


69 


177 5 


177 8 


2 


13 


17 78 


2 4 


128 3 


178 6. 


20 


1 4 


17 86 


68 


8 002 


179 5} 


20 


1 5 


17 95 


1 5 


156 8 


1773 


20 


1 6 


1773 


35 


0001 


177 9 


20 


1 7 


1779 


44 


178 3 


178 4 


20 


1 8 


1783 


4 6 


193 7 


193 0. 


20 


19 


17 84 


30 


00 5 


1897 


20 


20 


18 97 


69 


1901 


175 4 1 


20 


21 


17 54 


22 


126 


176 3, 


20 


2 2 


1763 


46 


176 6 


1918 


20 


2 3 


1918 


1 6 


1921 


192 6] 


20 


2 4 


1926 


4 6 


192 9 


1930 


20 


2 5 


1930 


60 


129 7 


12 6 4J 


20 


2 6 


19 29 


6 5 


126 5 


17 6 91 


20 


2 7 


17 69 


69 


1922 


1777. 


20 


2 8 


17 66 


60 


1919 


1923; 


2 


2 9 


1923 


64 


126 5 


177 ( 


20 


30 


1770 


69 


192 4 


177 7. 


20 


31 


17 77 


24 


128 5 


17 8 8' 


20 


3 2 


17 88 


3 5 


001 


1260 


20 


3 3 


1260 


30 


[49-b] 
8 002 


17 8 1. 


20 


3 4 


17 81 


67 


1789 


2 


3 5 


17 89 


20 


129 3 


179 6| 


20 


3 6 


1796 


1 5 


179 9 


175 2, 


20 


3 7 


17 52 


4 4 


175 5 


175 6. 


20 


3 8 


17 55 


60 


175 9 


1913 


20 


3 9 


1913 


1 9 


129 3 


192 7 


20 


4 


1927 


3 5 


0001 


193 3 


20 


41 


19 33 


10 


166 5 


192 


20 


42 


19 20 


21 


127 7 


178 


20 


4 3 


1780 


60 


129 3 


179 7 


20 


44 


1797 


11 


175 9 


176 4 


20 


4 5 


1764 


64 


127 7 


190 2 


20 


4 6 


1902 


69 


190 5 


190 9. 


20 


47 


17 56 


65 


129 3 


114 7 


20 


4 8 


1147 


3 5 


01 


175 3 


20 


4 9 


17 53 


69 


190 6 


190 9 


20 


50 


1909 


24 


126 7 


157 0) 


20 


51 


1570 


20 


127 8 


193 1A 


20 


5 2 


1931 


60 


8 2 


193 9 


20 


53 


19 39 


1 9 


8 001 


194 5 


20 


5 4 


19 45 


21 


125 


190 3 


20 


5 5 


1903 


61 


8 00 3 


1911 


20 


5 6 


1911 


1 9 


177 6 


194 


20 


57 


19 40 


60 


8 3 


114 9 


20 


5 8 


1149 


10 


16 5 6 


17 6 2 



Get B-B] ,bj store B-\ . 



Prepare to use the formula 
arctan (-x) ■ -arctan x 
at the end if B < 0. 



Test 49-b j go to special 
routines if |49-b| > 10 



Prepare for shift of 49-b. 

If b > 49, go to arctan 1/x. 
If b <; 46, set arctan x » x, 
go to end. 

If 46 < b £ 49, go to 
series evaluation 

If b > 49, calculate l/Bi . 



Shift to get x«B or x-l/B 
into fixed decimal form 

If jx| > a 28„ go to calc. 

V —V 



z ■ 



1 + 



-i i y».6 
xjy 



*]=£ 



1 + My* 



z m |xj. Prepare to 
add or arctan y. 



-xxi- 



20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 

20 
20 
20 
20 
20 
20 
20 
20 
20 
20 
20 



5 
6 
6 
6 
6 
6 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
7 
7 
7 
7 
7 



80 
81 



88 
8 9 
90 



00 
01 
2 
03 



8 
9 



1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



7 

7 

7 

9 

9 

1 

7 

78 

14 

90 

90 

91 

78 

791 

9 32 

267 

285 

9 43 

910 

934 

9 4-2 

1 5 

9 

1 

3 

4 

6 

8 64 
283 
2 64 
936 
937 



3 
6 

5 
4 
5 



58 
08 
28 
74 
75 
68 
01 
21 
22 
1 9 
24 
99 
59 
65 
05 



906 



8 5 
31 
80 
01 
18 
94 



60 
10 
1 9 
60 
10 

u 



60 
36 
11 
10 

U 

U 
65 



65 



8 00 
1 51 
8 00 

179 
1 70 



128 3 



8 



00 3 
00 3 



00 3 
00 3 
100 
00 3 
000 



731 

580 



15 7 07 9 



arc tan z ■ 

L 3 5 11 J 
(coefficients from LOG) 



I 



(correction at the end) 

arctan |x [-arctan z+£0 or arctan y] 

If b>49, arctan |x| -n/2- arctan l/|x|. 



If arctan II | > 1, round 
in yth decimal place, if 
< l f round in 8th, 



Normalize arc tan |5| 



C - arc tan B in radians 

For degrees, go to MFY 

by 360/2rr. 

For b > 59, arc tan |b| - n/2. 



144 5 

193 6 

194 1 
12 6 4 
1264 

178 1 

1910 

194 3 



128 5 
128 5 



r Constants 



4 41 9 5 



801 



57 
50 



2 95 7 



632 7 
005 
805 1 



J 



20 114 1143 19 1801 1267 Correction in arc tan z. 
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